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ABSTRACT 

We present a new, i^~selected, optical-to-near infrared photometric catalog of the Extended Chan- 
dra Deep Field South (ECDFS), making it publicly available to the astronomical community^ The 
dataset is founded on publicly available imaging, supplemented by original z' JK imaging data col- 
lected as part of the MUltiwavelength Survey by Yale-Chile (MUSYC). The final photometric catalog 
consists of photometry derived from UU^%BV RI z' J K imaging covering the full ^ y. ^ oi the 
ECDFS, plus H band photometry for approximately 80 % of the field. The ha flux limit for point- 
sources is k\^^^ — 22.0. This is also the nominal completeness and reliability limit of the catalog: 
the empirical completeness for 21.75 < K < 22.00 is > 85 %. We have verified the quality of the 
catalog through both internal consistency checks, and comparisons to other existing and publicly 
available catalogs. As well as the photometric catalog, we also present catalogs of photometric red- 
shifts and restframe photometry derived from the ten band photometry. We have collected robust 
spectroscopic redshift determinations from published sources for 1966 galaxies in the catalog. Based 
on these sources, we have achieved a {la) photometric redshift accuracy of A2;/(l + z) — 0.036, with 
an outlier fraction of 7.8 %. Most of these outliers are X-ray sources. Finally, we describe and release 
a utility for interpolating restframe photometry from observed SEDs, dubbed InterResfl. Particularly 
in concert with the wealth of already publicly available data in the ECDFS, this new MUSYC catalog 
provides an excellent resource for studying the changing properties of the massive galaxy population 
at z < 2. 

Subject headings: Catalogs — Techniques: Photometric — Galaxies: Observations — Galaxies: Distances 
and Redshifts — Galaxies: High-Redshift — Galaxies: Fundamental Parameters 



1. INTRODUCTION 

Over the past decade, multi-band deep-field imaging 
surveys have provided new opportunities to directly ob- 
serve the changing properties of the general, field galaxy 
population with lookback time. These new data, quanti- 
fying the star formation, stellar mass, and morphological 
evolution among galaxies, have led to new and funda- 
mental insights into the physical processes that govern 
the formation and evolution of galaxies. These advances 
have been made possible not only by the advent of a new 
generation of space-based and 8 m class telescopes, but 

^Imaging and spectroscopy data and catalogs are freely avail- 
able through the MUSYC Publi c Data Release webpage: 
|http://www.astro.yale.edu/MUSYC/j 

''InterRest can be do wnloaded from 

[W;tp://www. strw.leidenuniv.nl/~ent/InterRest Documenta- 
tion, including a complete walkthrough, is available from the same 
address. 



also the maturation of techniques for estimating redshifts 
and intrinsic properties like stellar masses from observed 
SEDs. These two developments have made it possible 
not only to go deeper — pushing to higher redshifts and 
probing further down the luminosity function — but also 
to consider many more galaxies per unit observing time. 
This has made possible the construction of large, repre- 
sentative, and statistically significant samples of galaxies 

spanning a large proportion of cosmic tirne; 

Th e Chandra Deep Field South (CDFS: GTacconi et all 
l2003 ) is one of the premier sites for deep field cosmologi- 
cal surveys (see Figure [1]). It is one of the most intensely 
studied region of the sky, with observations stretching 
from the X-ray to the radio, including ultraviolet, optical, 
infrared, and submillimeter imaging, from space-based as 
well as the largest terrestrial observatories. It has also 
become traditional for surveys targeting the CDFS, to 
make their data publicly available. As a direct result of 
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Fig. 1. — MUSYC in tlie ECDFS. — The greyscale image shows the new K band data. The solid black contour shows the area with useful 
photometry in all of UU^^BV RIz' JK in the MUSYC catalog. (Areas badly affected by bright stars in the z' band have been masked.) 
The catalog also includes H photometry for ~ 80 % of the field (s olid grey contour). For comparison, we als o show the area cove red by 
several other important (E)CDFS surveys: GEMS < Rix et al. ll2004l. dotted lines), the original Chandra CDFS (I Giacconi eran i2002!. short- 
dashed circle), the GOODS (Dickinson et al. 2002) HST ACS optical (Ught long-dashed rectangle) and ISAAC NIR (short-dashed region) 
imaging, the K20 survey (Cimatti et al. 2002, heavy long-dashed rectangle), and the HUDF (Beckwith et al. 2006, grey solid diamond). 
The FIREWORKS catalog (iWuvts et al. 2008) combines the GOODS ACS and ISAAC data with the UU^gBRVIz' data described in this 
paper for the central GOODS ISAAC region. SIMPLE (Damen et al., in prep.) will add very deep S pitzer IRAC imaging to the whole 
region shown here. A medium band NIR survey is also underway using the NEWFIRM instrument l iVan Dokkum et al.ll2059l) . At right, 
we show a detail of the K20 survey area (below), and futher detail of an approximately 2" X 2" area (above). 



this commitment to collaboration within the astronomi- 
cal community, the wealth of data available — in terms 
of both volume and quality — provide an exceptional 
opportunity to quantify the evolution of the galaxy pop- 
ulation out to high redshift. 

With this goal in mind, the key to gaining access to 
the z > 1 universe is near infrared (NIR) data. Most 
of the broad spectral features (e.g. the Balmer and 4000 
A breaks) on which modern SED-fitting algorithms rely 
are in the restframe optical; for z > 1, these features 
are redshifted beyond the observer's optical window and 
into the NIR. For this reason, we have combined exist- 
ing imaging of the Extended Chandra Deep Field South 
(ECDFS; see Figure [T]) with new optical and NIR data 
taken as part of the MUltiwavelength Survey by Yale- 
Chile (MUSYC). 

The primary objective of MUSYC is to obtain deep op- 
tical imaging and spectroscopy of four ^ x ^ South- 
ern fields, providing parent catalogs for followup with 
ALMA. Coupled w ith the optical {UBVRIz') imaging 
program (jCawiser et al. 2006a) , there are two NIR com- 
ponents t o the MUSYC pro ject: a deep component 
[K < 23.5: lQuadri et al.ll2007t l. targeting four 10 x 10 □" 
regions w ithin the MUSYC fields, and a wide component 
{K < 22: lBlanc et~alll2008l this work) covering three of 



the i X i MUSYC fields in their entirety. These 
data are intended to allow, for example, the restframe- 
UV selection of gal axies at z > 3 us ing the Lyman 
break technique fe.g. ISteidel et al1ll996[ ). the restframe- 
optical selection of galaxies at z > 2 usin g the Distant 
Red Galaxy (DRG) criterion (jFranx et al.|[200 3). and the 
color -selection of z ^,1 -4 galaxies using the BzK crite- 
rion (|Daddi et al.l[200l . 

In the ECDFS, the broadband imaging data have 
been supplemented by a narrow-band imaging s urvey, 
targeting Ly-a emitters at z — 3.1 (Gawiser et al.l 
|2O06b; Gronwall et al. 2007), and a spectroscopic sur- 
vey (jTreister et al.. .2008.) targeting Xray sources froni 
the 250 ks ECDF S Xray catalog ()Lehmer et al.l 120051 
iVirani et al.|[2006[ ). Further, the Spitzer IRAC/MUSYC 
Public Legacy in the ECDFS (SIMPLE; M Damen et 
al., in prep.) project has obtained very deep IRAC 
imaging across the full ECDFS. There is also a deep 
medium band optical survey underway (Cardamone et 
al., in prep.), and a plann ed medium band NIR survey 
()Van Dokkum et al.|[2009( ). 

This paper describes the MUSYC wide NIR-selected 
catalog of the ECDFS (which we will from now on refer to 
as 'the' MUSYC ECDFS catalog, despite the existence of 
several separate MUSYC catalogs, as described above), 



A isT-Selected Catalog of the ECDFS from MUSYC 



3 



TABLE 1 

Summary of the data comprising the MUSYC ECDFS catalog 
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B 


4605 
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-0.12 


19.29 


1012 


1'.'03 


26.9 


15153 


0.852 


8223 


8322 


880 


V 


5383 


895 


-0.01 


29.06 


1022 


0"95 


26.6 


15154 


0.863 


8370 


8463 


891 


R 


6520 


1600 


+0.19 


24.35 


1017 


0"88 


26.3 


15148 


0.894 


8647 


8758 


897 


I 


8642 


1500 


+0.51 


9.60 


977 


0^98 


24.8 


15128 


0.826 


8456 


8545 


897 


z' 


9035 


995 


+0.54 


1.30 


996 


1('13 


24.0 


13972 


0.751 


8043 


8000 


897 


J 


12461 


1620 


+0.93 


1.33 


906 


< 1'.'49 


23.1 


14580 


0.683 


7894 


7859 


896 


H 


16534 


2960 


+1.40 


1.00 


560 


< 1'.'22 


23.1 


10518 


0.579 


7005 


6313 


692 


K 


21323 


3310 


+1.83 


1.00 


906 


< 1'.'05 


22.4 


14355 


0.695 


8782 


8911 


897 



Note. — For each band (Col. 1) that has gone into the MUSYC ECDFS catalog, we give the effective wavelength (Col. 2), the filter 
FWHM (Col. 3), and the apparent magnitude of Vega, in the AB system (i.e. the conversion factor between the AB and Vega magnitude 
systems, Col. 4). We also give the mean integration time (Col. 5) for each image, the effective imaging area (defined as the region receiving 
more than 75 % of the nominal integration time, Col. 6), and the final effective seeing (FWHM, Co l. 7). The 5(t limiting depths given in 
Col. (8) are as measured in 2'.'5 diameter apertures on the 1'.'5 FWHM PSF-matched images (see i]3.2l for a point source, these can be 
translated to total magnitudes by subtracting 0.45 mag. Note that, whereas the optical data are taken in single pointings, the final NIR 
images are mosaics of many pointings. Note that the central ~ 10 X 10 of the field received an extra three hours' integration time in 
the H band; these data are approximately 0.3 mag deeper than the figure quoted above. Col. (9) gives the number of K detections that 
useful coverage {i.e. an effective weight, w, of 0.6 or greater) in each band ; Col. (10) gives the fraction of those objects that have > 5cr 
detections. Both of these columns refer to the full catalog. Col. (11) gives the number of objects in the main science sample (i^Tot < 22, 
K S: N > 5, wb > 0.6, w^i > 0.6, wk > 0.75) with > 5cr detections; Col.s (12) and (13) give the numbers of stars and galaxies separately 
(see EIJ. 



and makes it publicly available to the astronomical com- 
munity. A primary scientific goal of the wide NIR com- 
ponent of the survey is to obtain statistically significant 
sample s of massive galax ies at z < 2. In a companion 
paper (jTaylor et al. 2009, hereafter Paper II), we will use 
this dataset to quantify the z <2 color and number den- 
sity evolution of massive galaxies in general, and in the 
relative number of red sequence galaxies in particular. 

The MUSYC ECDFS dataset is founded on existing 
and publicly available imaging, supplemented by original 
optical (z') and NIR (JK) imaging. Apart from the JK 
imaging, all these data have been described elsewhere. 
Accordingly, the data reduction and calibration of the 
new JK imaging is a prime focus of this paper. How- 
ever, when it comes to constructing panchromatic cata- 
logs with legacy value from existing datasets, the whole 
is truly more than the sum of parts: ensuring both abso- 
lute and relative calibration accuracy is paramount. We 
have invested substantial time and effort into checking 
all aspects of our data and catalog, using both simulated 
datasets, and through comparison to some of the many 
other existing (E)CDFS catalogs. 

The structure of this paper is as follows: we de- 
scribe the acquisition and basic reduction of the MUSYC 
ECDFS broadband imaging dataset in !j2j The pro- 
cesses used to combine these data into a mutually con- 
sistent whole are described in fJSl In SJH we describe the 
construction of the photometric catalog itself, including 
checks on the completeness and reliability, and on our 
ability to recover total fluxes. We present external checks 
on the astrometric and photometric calibration in [JS] Af- 
ter a simple comparison of our catalog to other NIR- 
selected catalogs in !|6l we describe our basic analysis of 
the multi-band photometry in fjTl including star/galaxy 
separation, and the derivation of photometric redshifts, 
as well as the tests we have performed to validate our 
analysis. In fJHl we introduce InterRest; a new utility for 



interpolating restframe fiuxes. This utility is also being 
made public. Additionally, in Appendix [Aj we describe 
a compilation of 2213 robust spectroscopic redshift de- 
terminations for objects in the MUSYC ECDFS catalog. 

Throughout this work, all magnitudes are expressed 
in the AB system; the only exception to this is §5.21 
where it will be convenient to adopt the Vega sys- 
tem. Where necessary, we assume the concordance cos- 
mology; viz. f^m = 0.3, r^A — 0.7, Hq — 1.0, and 
Hq — 70 km s^^ Mpc^^. When discussing photometric 
redshifts, we will characterise random errors in terms of 
the NMAD^ of Az/(1 + z); we will abbreviate this quan- 
tity using the symbol a^. 

2. DATA 

This section describes the acquisition of the imag- 
ing data comprising the MUSYC ECDFS dataset; the 
vital statistics of these data are given in Table [TJ 
Of these data, only the z' JK are original; the WFI 
UU^fiBV RI imaging has been reduced and described by 
Hildebrandt ct al. (2006), and the SofI H band data by 
iMov et al.l (l2003l) Further, the origin al z' data have been 
reduced as per lGawiser et al.l (|2006al ) for the MUSYC op- 
tical {BVR — selected) catalog. We have therefore split 
this section between a summary of the data that are de- 
scribed elsewhere f §2.1|) . and a description of the new 
ISPI JK imaging (i i2.2|) . Note that what we refer to as 
the K band is really a 'if short' filter; we have dropped 
the subscript for convenience. For a complete description 
of the other datasets, the reader is referred to the works 
cited above. 

2.1. Previously Described Data 

^ Here, NMAD is an abbreviation for the Normalized Median 
Absolute Deviation, and is defined as 1.48 X raed[x — med(x)]; the 
normalization factor of 1.48 ensures that the NMAD of a Gaussian 
distribution is equal to its standard deviation. 
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2.1.1. The WFI Data — UUsgBVRI Imaging from the 
ESO Archive 

iHildebrandt et al] (|2006) have collected all (up un- 
til December 2005) archival UUssBRVI^ imaging data 
taken using the Wide F ield Imager (WFI, 0'.'238 pix'^; 
iBaade et al.lfl998l fl999l) on the ESO MPG 2.2 m tele- 
scope for the four fi elds that make up t he ESO Deep 
Public Survey (DPS: Arnouts e t al. 2001). In addition 
the original DPS ECDFS data (DPS field 2c), this com- 
bined dataset includes WFI commi ssioning da. t a, th e 
data from the COMBO- 17 survey (|Wolf et al.l [200I . 
and observations f r om se ven other observing programs. 
IHildebrandt et ahl (|2006f ) have pooled and re-reduced 
these data using the automa ted THELI pipeline de- 
scribed by lErben et al.l ()2005[ ) under the moniker Ga- 
BoDS (Garching Bonn Deep Survey). The final prod- 
ucts are publicly available through the ESO Science 
Archive Facility.^ The final ima ge quality of these images 
is 0'.'9— T.'l FWHM. .Hildebrandt et all (j2006) estimate 
that their basic calibration is accurate to better than 
^ 0.05 mag in absolute terms, and that, based on color- 
color diagrams for stars, the relative or cross-calibration 
between bands is accurate to < 0.1 mag for all images. 

2.1.2. The Mosaic-II data — Original z Imaging 

We have supplemented the WFI optical data with 
original z' band imaging taken usi ng Mosaic-II cam- 
era (0'.'267 pix-^ iMuller et al.lll998i ) on the CTIO 4m 
Blanco telescope. The data acquisition strategy is the 
same as for the opti cal data in other MUSYC fields 
(|Gawiser et al.ll2006aD : the ECDFS data were taken in 
January 2005. The final integration time was 78 min- 
utes, with an effective seeing of I'.'l FWHM, although 
we note that the PSF does have broad, non-Gaussian 
'wings'. The estimated u ncertainty in the pho tometric 
calibration is < 0.03 mag (|Gawiser eirall[2006ah . 

2.1.3. The SofI Data — H Imaging Supporting the ESO 

DPS 

We include the H band data described by iMov et ahl 
(HOOI), which was taken to complement t he original DPS 
WFI optical data and SofI NIR data (jVandame et al.l 
[2OOII: IOlsCT"eraI][2006l) . This dataset covers approxi- 
mately 80 % of the ECDFS, consisting of 32 separate 
4'.'9 X 4'.'9 pointings, and were obtaining using SofI (0'.'288 
pix-i; iMoorwood et al1ll998f ) on the ESO NTT 3.6 m 
telescope. The data were taken as a series of dithered 
(or 'jittered') 1 min exposures, totaling 60 min per point- 
ing; the central four fields received an extra 3 hours' ex- 
posure time. We received these data (P aulin e Barmby, 
priv. comm.) reduced as described bv lMov et al.l (20031 ): 
i.e., as 32 separate, unmosaicked fields. The effective 
seeing in each po inting varies from 0'.'4 to 0'.'8 FWHM. 
iMov et all ()2003[ ) found that their photometric zeropoint 

^ Two separate WFI U filters have been used. The first, 
ESO#877, which we refer to as the U filter, is slightly broader 
than a Broadhurst U filter. This filter is known to have a red leak 
beyond 8000 A. The second filter, ESO#841, which we refer to as 
C/38, is something like a narrow Johnson U filter. There is, unfor- 
tunately no clear convention for how to refer to these filters; for 
instance. lArnouts et al.l ||200H ) refer to what we call the U and U3S 
as U' and U, respectively. 

^ http:/ /archive. eso. org/cms/eso-data/data-packages/gabods- 
data-release- version- 1.1-1/ 



solution varied by < 0.04 mag over the course of a night; 
they offer this as an upper limit on possible calibration 
errors. Further, in comparis on to the Los Campanas In- 
frared Survey (LCIRS; Che n et al.l[200ll . and the vO.5 
(April 200 2) release of the GOODS ISAAC photometry, 
iMov et all (,200 3) found their calibration to be 0.065 mag 
brighter, and 0.014 mag fainter, respectively. 

2.2. The ISPI Data — Original JK Imaging 

The new MUSYC NIR imaging consists of two mo- 
saics in the J and K bands, each made up of 3 x 3 
pointings, and covering approximately 950 □'. The data 
were obtained using the Infrared Sideport Imager (ISPI 
- Probst et al. 2003; Van der Bhek et al. 2004) on the 
CTIO Blanco 4m telescope. ISPI uses a 2048 x 2048 
pix HgCdTe HAWAII-2 detector, which covers approxi- 
mately 10'.'5 x 10'.'5 at a resolution of « 0'.'3 pix'^ The 
aim was to obtain uniform J and K coverage of the full 
i X i of the ECDFS to - 80 minutes and ~ 60 min- 
utes, respectively; our target (5ct, point source) limiting 
magnitudes were J « 22.5 and K k, 22. 

The data were taken over the course of fifteen nights, 
in four separate observing runs between January 2003 
and February 2004. In order to account for the bright 
and variable NIR sky (~ 10000 times brighter than 
a typical astronomical source of interest, varying on 
many-minute timescales), the data were taken as a 
series of short, dithered exposures. A non-regular, 
semi-random dither pattern within a 45" box was used 
for all but three sub-fields; these three earliest pointings 
were dithered in regular, 10" steps. An exposure of 
4 X 15s {i.e., 4 individual integrations of 15 seconds, 
coadded) was taken at each dither position in K; in J, 
exposures were typically 1 x 100s. 

Conditions varied considerably over the observing cam- 
paign, with seeing ranging from < 0'.'7 to > 1"5 FWHM. 
All nine K band pointings were observed under good con- 
ditions ( < I'.'O FWHM). However, observing condititions 
were particularly bad for two of the nine J pointings; the 
final effective seeing of both the South and Southwest 
pointings are nearer to 1'.'5 FWHM. 

For each of the subfields comprising the MUSYC ISPI 
coverage of the ECDFS, the data reduction pipeline is 
essentially the sa me as for the o t her M USY C NIR imag- 
ing, d escribed by lOuadri et al.l (|2007D and iBlanc etahl 
120081 ) , foll o wing the same basic strategy as, e.g., 
Labbi Tet al.l (|2003f ). The data reduction itself was per- 
formed using a modified version of the IRAF package 
xdimsum."' 

2.2.1. Dark Current and Flat Field Correction 

The ISPI detector has a non- negligible dark current. 
To account for this, nightly 'dark flats' were constructed 
by mean combining (typically) ten to twenty dark expo- 
sures with the appropriate exposure times; these 'dark 
flats' are then subtracted from each science exposure. 
These dark flats show consistent structure from night to 

IRAF is distributed by the National Optical Astron- 
omy Observatories, which are operated by the Associ- 
ation of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science 
Foundation. The xdimsum package is available from 

http: / / iraf.noao.edu /iraf/ ftp /iraf / extern- v2 12/ xdimsum020806 . 
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night, but vary somewhat in their actual levels. Note 
that this correction is d one before fiat-fie lding and/or 
sky subtraction (see also lBlanc et al.ll2008[ ). 

Flat field and gain/bias corrections {i.e., spatial varia- 
tions in detector sensitivity due to detector response, op- 
tic throughput, etc.) were done using dome-flats, which 
were constructed either nightly or bi-nightly. These flats 
were constructed by taking a number of exposures with 
or without a lamp lighting the dome screen. Each flat- 
field was constructed using approximately ten 'lamp on' 
and 'lamp off' exposures, mean combined. In order to 
remove background emission from the 'lamp on' image, 
we subtract away the 'lamp off' image, to leave only 
the light reflected from off of the dome screen (see also 
iQuadri et al.ll2007l ). These flats are very stable night to 
night, with some variation between different observing 
runs. 

2.2.2. Sky Subtraction and Image Combination 

Because the NIR sky is bright, non-uniform, and vari- 
able, a separate sky or background image must be sub- 
tracted from each individual science exposure. The basic 
xdimsum package does this in a two-pass procedure. In 
the first pass, a background map is constructed for each 
individual science image by median combining a sequence 
of (typically) eight dithered but temporally continguous 
science exposures: typically the four science images taken 
immediately before and after the image in question. In 
the construction of this background image, a 'sigma clip- 
ping' algorithm is used to identify cosmic rays and/or 
bad pixels, which are then masked out. The resultant 
background image (which at this stage may be biased 
by the presence of any astronomical sources) is then sub- 
tracted from the science image to leave only astronomical 
signal. The sky subtracted images are then shifted to a 
common reference frame using the positions of stars to 
refine the geometric solution, undoing the dither, and 
then mean combined, again masking bad pixels/cosmic 
rays. This combined image is used to identify astro- 
nomical sources, using a simple thresholding algorithm. 
This process is repeated in the second 'mask pass', with 
the difference that astronomical sources are now masked 
when the b ackground map is co nstructed. 

Following IQuadri et all (|2007t l , we have made several 
modifications to the basic xdimsum algorithm in order to 
improve the final image quality. We have constructed an 
initial bad pixel mask using the flat-field images. Fur- 
ther, each individual science exposure is inspected by 
eye, and any 'problem' exposures (especially those show- 
ing telescope tracking problems or bad background sub- 
traction) are discarded; artifacts such as satellite trails 
and reflected light from bright stars are masked by hand. 
These masks are used in both the first pass and mask 
pass. 

Persistence is a problem for the ISPI detector: as 
a product of detector memory, 'echoes' of particularly 
bright objects linger for up to eight exposures. For this 
reason, we have also modified xdimsum to create sepa- 
rate masks for such artifacts; these masks are used in 
the mask pass. Note that for the three subfields (includ- 
ing the Eastern K pointing) observed using a regular, 
stepped dither pattern, this leads to holes in the cover- 
age near bright objects: the 'echoes' fall repeatedly at 
certain positions relative to the source, corresponding to 



the regular steps of the dither pattern. At worst, cover- 
age in these holes is ~ 25% of the nominal value. 

Even after sky-subtraction, large-scale variations in the 
background were apparent; these patterns were different 
and distinct for each of the four quadrants of the images, 
corresponding to ISPI's four amplifiers. To remove these 
patterns, we have fit a 5th order Legendre polynomial 
to each quadrant separately, using 'sigma clipping' to 
reduce the contribution of astronomic al sources, and the n 
simply subtracted this away (see also lBlanc et al.|[2008[ ). 
This subtraction is done immediately after xdimsum's 
normal sky-subtraction. 

In the final image combination stage, we adopt 
a weighting scheme designed to optimize signal-to- 
noise for point sou rces (see, e.g.. iGawiser et al]l2006al : 
IQuadri et al.ll2007f) . At the end of this process, xdim- 
sum outputs a combined science image. Additionally, 
xdimsum outputs an exposure or weight map, and a map 
of the RMS in coadded pixels. Note that although this 
RMS map is not accurate in an absolute sense, it does 
do an adequate job of mapping the spatial variation in 
the noise; see i )4.6l 

2.2.3. Additional Background Subtraction 

The sky subtraction done by xdimsum is imperfect; a 
number of large scale optical artifacts (particularly re- 
flections from bright stars and 'holes' around very bright 
objects) remain in the K images as output by xdimsum. 
Using these images, in the object detection/extraction 
phase, we were unable to find a combination of SExtrac- 
tor background estimation parameters (viz. BACK_SIZE 
and BACK_FILTERSIZE) that was fine enough to map 
these and other variations in the background but still 
coarse enough to avoid being influenced by the biggest 
and brightest sources. This led to significant incomplete- 
ness where the background was low, and many spurious 
sources where it was high. We were therefore forced to 
perform our own background subtraction, above and be- 
yond that done by xdimsum. 

This basic idea was to use SExtractor 'segmentation 
maps' associated with the optical (BVR^) and NIR (K) 
detection images to mask real sources. In particular, 
the much deeper BVR stack includes many faint sources 
lying below the K detection limit. To avoid the con- 
tributions of low surface brightness galaxy 'wings', we 
convolved the combined [BV R+K) segmentation maps 
with a 15 pix (4") boxcar filter to generate a 'clear 
sky' mask. Using this mask to block flux from astro- 
nomical sources, we convolved the science image with 
a 100 pix (26'.'7) FWHM Gaussian kernel to generate a 
new background map; this was then subtracted from the 
xdimsum-generated science image. 

Note that the background subtraction discussed above 
is important only in terms of object detection; back- 
ground subtraction for photometry is discussed in tj4.3l 
While this additional background subtraction step re- 
sults in a considerably flatter background across the de- 
tection image, it does not significantly or systematically 
alter the measured fluxes of most individual sources. 

2.2.4. Photometric Calibration 

^ Here, by BVR, we are referr ing to the combined B + V + R op- 
tical stack used for detection bv lGawiser et al.l l|2006bl l in the con- 
struction of the MUSYC optically-selected catalog of the ECDFS. 
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Because not all pointings were observed under photo- 
metric conditions, we have secondarily calibrated each 
NIR pointing separately with reference to the 2MASS 
(jCutri et all l2003l : ISkrutskie et~an l2006f) Point Source 
Catalog.^ Taking steps to exclude saturated, crowded, 
and extended sources, we matched ISPI magnitudes mea- 
sured in 16" diameter apertures to the 2MASS catalog 
'default' magnitude (a 4" aperture flux, corrected to to- 
tal assuming a point-source profile). For each subfield, 
the formal errors on these zeropoint determinations are 
at the level of 1 — 2 percent. The uncertainty is domi- 
nated by the 2MASS measurement errors, and are high- 
est for the central pointing where there are only 6 — 8 
useful 2MASS~detected point sources. For comparison, 
the formal 2MASS estimates for the level of systematic 
calibration errors is < 0.02 mag. 

3. DATA COMBINATION AND 
CROSS-CALIBRATION 

This section is devoted to the combination and cross- 
calibration of the distinct datasets described in the pre- 
vious section into a mutually consistent whole. In 
we describe the astrometric cross-calibration of each of 
the ten images, including the mosaicking of the NIR 
data. We describe and validate our procedure for PSF- 
matching each band in ij3.2l 

3.1. Astrometric Calibration and Mosaicking 

To facilitate multi-band photometry, each of the final 
science images is transformed to a common astrometric 
reference frame: a North-up tangential plane projection, 
with a scale of 0'.'267 pix^^. This chosen reference frame 
corresponds to the stacked BVR image used as the de- 
tection imag e for the optically-sele cted MUSYC ECDFS 
catalog (see iGawiser et al1l2006al lbl). based on an early 
reduction of the WFI data. 

Whereas WFI and Mosaic-II are both able to cover the 
entire ECDFS in a single pointing, the Sofl and ISPI cov- 
erage consists of 32 and nine subfields, respectively. For 
these bands, each individual subfield was astrometrically 
matched to the BVR reference image using standard 
IRAF/PyRAF tasks. For the ISPI data, each subfield 
is then combined, weighted by S:N on a per pixel ba- 
sis, in order to create the final mosaicked science image. 
(Note that individual subfields are also 'PSF-matched' 
before mosaicking - see ^13.21 ) 

One severe complication in this process is that expo- 
sure/weight maps were not available for the Sofl imag- 
ing. We have worked around this problem by construct- 
ing mock exposure maps based on estimates of the per 
pixel RMS in each science image. Specifically, we calcu- 
late the biweight scatter in rows and columns: <7b(x) and 
(Tb(j/). The effective weight for the pixel (x, y) is then 
estimated as [cB{x)(JB{y)]^'^ ■ The map for each subfield 
is normalized so that the median weight is 1 for those 
pointings that received 1 hour's integration, and 4 for 
the four cent ral pointings. 

In line with'Qua dri et all ()2007[ ). we found it necessary 
to fit a high order surface (wz., a 6th order Legendre 
polynomial, including x and y cross terms) to account for 
the distortions in the ISPI focal plane. For the Sofl data, 

^ Available electronically via GATOR: 

|http://irsa.ipac. caltech.edu/applications/Gator/ 1 
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Fig. 2. — The astrometric registration of the (from top to bot- 
tom) Iz'JH images (obtaining using WFI, Mosaic-II, ISPI, and 
Soft, respectively), relative to the K detection image — In each 
panel, vectors give the biweight mean positional offset between 
the two images in 2.'5 X 2.'5 cells, based on all K < 22 sources; 
the greyscale gives the biweight variance. Systematic astrometric 
shears in individual images images are typically much less than a 
pixel. 
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a 2nd order surface was sufficient, although we did find 
it necessary t o revis e the initial astrometric calibration 
bv lMov et all (|200l . 

As an indication of the relative astrometric accuracy 
across the whole dataset, Figure [2] illustrates the differ- 
ence between the positions of all K < 22 sources mea- 
sured from the K band, and those measured in each of the 
Rz'JH bands (observed using, in order, WFI, Mosaic-II, 
ISPI, and SofI). Systematic 'shears' between bands are 
typically much less than a pixel. Comparing positions 
measured from the registered R and K band images, av- 
eraged across the entire field, the mean positional offset 
is 0'.'15 (0.56 pix). Looking only at the x/y offsets, we 
find the biweight mean and variance to be O'.'OS (0.11 pix) 
and O'.'S (1.1 pix), respectively. 

3.2. PSF Matching 

The basic challenge of multi-band photometry is ac- 
counting for different seeing in different bands, in order 
to ensure that the same fraction of light is counted in each 
band for each object. We have done this by matching the 
PSFs in each separate pointing to that with the broad- 
est PSF. Of all images, the South- Western J pointing 
has the broadest PSF: I'.'S FWHM. This sets the limiting 
seeing for the multiband SED photometry. Among the 
K pointings, however, the worst seeing is I'.'O FWHM; 
this sets the limiting seeing for object detection, and the 
measurement of total K magnitudes (see ij4.1l and ij4.3p . 
We have therefore created eleven separate science im- 
ages: one I'.'S FWHM image for each of the ten bands to 
use for SED photometry, plus a I'.'O FWHM K image for 
object detection and the measurement of total K fluxes. 

The PSF-matching procedure is as follows: for each 
pointing, we take a list of SED-classified stars from the 
COMBO-17 catalog; these objects are then used to con- 
struct an empirical model of the PSF in that image, using 
an iterative scheme to discard low signal-to-noise, ex- 
tended, or confused sources. Our results do not change 
if we begin with Bz'K selected stars, or GEMS point 
sources. We then use the IRAF/PyRAF task lucy (an 
implementation of the Lucy-Richardson deconvolution 
algorithm, and part of the STSDAS package^) to deter- 
mine the convolution kernel required to 'degrade' each 
subfield to the limiting effective seeing. Finally, the con- 
volution is done using standard tasks. Note that each 
of the NIR subfields is treated individually, prior to mo- 
saicking. 

In order to quantify the random and systematic errors 
resulting from imperfect PSF matching. Figure [3] shows 
the relative difference between the curves of growth of in- 
dividual point sources across nine of our ten bands, after 
matching to the target 1'.'5 FWHM PSF. In this Figure, 
we compare the growth curves of many bright, unsat- 
urated, isolated point sources as a function of aperture 
diameter; specifically, we plot the relative difference be- 
tween the normalized growth curves in each band, com- 
pared to the median K band growth curve. Within each 
panel, the circles represent the median growth curve in 
each band (zero for the K band by construction), and 
the large and small error bars represent the 33/67 and 
5/95 percentiles, respectively. 

STSDAS is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 
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Fig. 3. — Relative deviations in the curves of growth for point- 
sources in each of nine bands, from four different instruments, af- 
ter PSF matching (I'.'S FWHM) — Each panel shows the relative 
differences between the normalized growth curves of bright, unsat- 
urated, isolated point sources, plotted as a function of aperture 
diameter. Circles show the median of all growth curves in each 
band; large and small error bars show the 33/67 and 5/95 per- 
centiles, respectively. The growth curves in different bands are all 
normalized with respect to the K band median; the systematic er- 
rors in the K panel are thus zero by construction. For our smallest 
apertures (2'' 5), systematic offsets due to imperfect PSF matching 
are at worst 0.006 mag; random errors, due to, for example, spatial 
variation of the PSF, are < 0.03 mag. 

After PSF matching, there are signs of spatial varia- 
tions in the FWHM of the J and K PSFs at the few per- 
cent level, particularly towards the edges of each point- 
ing. But since the scatter in these plots represents both 
real spatial deviations in the PSF, as well as normaliza- 
tion errors, these results can thus can be taken as an 
upper limit on the random PSF-related photometric er- 
rors. Looking at the z'-band panel, it is possible that the 
broad z' band PSF wings are important at the < 0.005 
mag level for 2'.'5 — 5'.'0. Note that the smallest apertures 
we use are 2'.' 5 in diameter — for these apertures, ran- 
dom errors due to imperfect PSF matching are typically 
< 0.03 mag, and systematic errors are at worst 0.006 
mag. 

4. DETECTION, COMPLETENESS, 
PHOTOMETRY, AND PHOTOMETRIC ERRORS 

In this section, we describe our scheme for building 
our multi-colour catalog of the ECDFS; a summary of 
the contents of the final photometric catalog i s given in 
Tabled We rely on SExtractor (Bcrti n fc Arn outs 1996|) 
for both source detection and photometry; in i|4.1l we de- 
scribe our use of SExtractor, and we quantify catalog 
completeness and reliability in t j4.2| There are two sep- 
arate components to the reported photometry for each 
object: the total K flux, which is discussed in t j4.3l and 
! j4.41 and the ten band SED, which is discussed in tj4.5l 
Finally, in ij4.61 we describe the process by which we have 
quantified the photometric measurement uncertainties. 
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TABLE 2 

Summary of the Contents of the Photometric Catalog 



Column No. Column Title 



Description 



1 

2, 3 
4 

5, 6 
7 

8—27 
28 

29, 30 



31, 32 
33, 34 
35—37 

38—47 
48 

49, 50 
51 

52—54 
55, 56 
57, 58 



id 

ra, dec 
field 
X, y 
ap_col 

U_colf, U_colfe, etc. 
ap_tot 

K.totf, K.totfe 



K_4arcsecf, K_4arcsecfe 
K_autof, K_autofe 
Kr50, Keps, Kposang 

Uw, etc. 
id_sex 

Ldeblcndl f_deblend2 
starjlag 

z_spec, qLspec, spec_class 
source, nsources 
qz_spec, spec_flag 



Object identifier, beginning from 1 

Right ascension and declination (J2000), expressed in decimal degrees 
An internal MUSYC field identifier (ECDFS=8) 
Center of light position, expressed in pixels 

Effective diameter {i.e. ^/4ttA, where A is the aperture area), in arcsec; we use the larger of 
SExtractor's ISO aperture and a 2''5 diameter aperture to measure colors (see i]4.5|l 
Observed flux," with the associated measurement uncertainty, in each of the 
UUssBV RI z' ,J H K bands, as measured in the 'color' aperture 

Effective diameter of the AUTO aperture, on which the total K flux measurement is based 
Total K flux — based on SExtractor's AUTO mea sure ment, with corrections applied for 
missed flux and background over-subtraction (see il4.3ll — and the associated measurement 
uncertainty, which accounts for correlated noise, random back ground subtraction errors, spa- 
tial variations in the noise, Poisson shot noise, etc. (see il4.6|l 

K flux, as measured in a 4" aperture, with the associated measurement uncertainty 
K flux within SExtractor's AUTO aperture, with the associated measurement uncertainty 
Morphological parameters from SExtractor, measured from the I'.'O FWHM K image; viz., the 
half-light radius (where the 'total' light here is the AUTO flux), ellipticity, and position angle 
Relative weight in each of the UUigBV RIz' JHK bands. 

The original SExtractor identifler,"^ for use with the SExtractor generated segmentation map 
Deblending flags from SExtractor, indicating whether an object has been deblended, and 
whether that olaject's photometry is significantly affected by a near neighbor, respectively 
A flag indicating whether an object's Bz' K color suggests its being a star (see t|7.1l l 
Spectroscopic redshift determination, if available, along with the associated quality flag and 
spectral classification, if given. 

A code indicating the source of the spectroscopic redshift, and the number of agreeing deter- 
minations 

A figure of merit, derived from the MUSYC photometry, for the spectroscopic redshift deter- 
mination (see Appendix |XJ , and a binary fiag indicating whether the spectroscopic redshift 
is considered 'secure' 



Note. — All fluxes are given in such a way that they can be transformed to AB magnitudes using a zeropoint of 25; in other words, 
fluxes are given in units of 0.363 mjy. For all but the z' and H bands, this is essentially the exposure time, normalized by the nominal 
values given in Table [T] For the H band, this value is derived from the mock exposure map described in il3.1l the z' band, this is a binary 
flag indicating whether the z' photometry is significantly affected by light from a nearby bright star. Recall that we have excised aU 
detections with an effective exposure time of less than 12 minutes from the catalog output by SExtractor. 



4.1. Detection 

Source detection and photometry for each band was 
performed using SExtractor in dual image mode; that 
is, using one image for detection, and then performing 
photometry on a second 'measurement' image. In all 
cases, the I'.'O FWHM K band mosaic (see §3.2p was used 
as the detection image; since flexible apertures are al- 
ways derived from the detection image, this assures that 
the same apertures are used for all measurements in all 
bands. 

As a standard part of the SExtractor algorithm, the 
detection image is convolved with a 'filter' function that 
approximates the PSF; we use a 4 pix (~ I'.'O) FWHM 
Gaussian filter. We adopt an absolute detection thresh- 
old equivalent to 23.50 mag / □" in the filtered image, re- 
quiring 5 or more contiguous pixels for a detection. Since 
we have performed our own background subtraction for 
the NIR images (see i j2.2.3p . we do not ask SExtractor to 
perform any additional background subtraction in the de- 
tection phase. For object deblending, we set the parame- 
ters DEBLEND_NTHRESH and DEBLEND.MINCONT 
to 64 and 0.001, respectively. These settings have been 
chosen by comparing the deblended segmentation map 
for the K detection image to the optical BVR detection 
stack, which has a considerably smaller PSF. 



Near the edges of the observed region, where cover- 
age is low, we get a large number of spurious sources. 
We have therefore gone through the catalog produced 
by SExtractor, and culled all objects where the K ef- 
fective weight, wk, is less than 0.2 (equivalent to < 12 
minutes per pointing). This makes the effective imaging 
area 953 □". Further, we find that a large number of 
spurious sources are detected where there are 'holes' in 
the coverage map (a product of the regular dither pat- 
tern used for the earliest eastern and northeastern tiles; 
see ^3.11 ) To avoid these spurious detections, for scien- 
tific analysis we will consider only those detections with 
an Wk > 0.75 (equivalent to ~ 45 minutes per pointing) 
or greater.^ This selection reduces the effective area of 
the catalog to 887 □". 

4.2. Completeness and Reliability 

In order to estimate the catalog completeness, we have 
added a very large number of simulated sources into the 
I'.'O FWHM detection image, and checked which are re- 

* In other words, the catalog is based on the area that received 
the equivalent of > 12 min integration, but our scientific analysis is 
based on those oljjects that received > 45 min integration. While 
objects with 0.2 < wk < 0.75 are given in the catalog, we do not 
include them in our main science sample, because of the poorer 
completeness and reliability among these objects. 
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covered by SExtractor, using the same settings as 'hve' 
detection. The completeness is then just the fraction 
of inputed sources which are recovered, as a function of 
source size and brightness. We adopted a de Vaucouleurs 
(i?^/'*-law) profile for all simulated sources, each with a 
half-light radius, ReS, between 0" {i.e. a point source) 
and 3", an ellipticity of 0.6, and total apparent K mag- 
nitude in the range 18 — 23 mag. We truncate each ob- 
ject's profile at 8 i?cff • No more than 750 artificial galax- 
ies were added at any given time, corresponding to 3-5 
% increase in the number of detected sources. Simulated 
sources were placed at least 13'.'35 (50 pix) away from 
any other detected or simulated source; these complete- 
ness estimates therefore do not account for confusion. 

The results of this exercise are shown in Figured which 
plots the completeness as a function of size and bright- 
ness. For point sources, we are 50 %, 90 %, and 95 % 
complete for K — 22.4, 22.2, and 22.1 mag, respectively. 
At a fixed total magnitude, the completeness drops for 
larger, low surface brightness objects. At K = 22, the 
nominal completeness limit of the catalog, we are in fact 
only 84 % complete for Rcs — 0'.'4, assuming an R^^^ 
profile. Note that we detect quite a few objects that 're- 
ally' lie below our formal (surface brightness) detection 
limit: just as noise troughs can 'hide' galaxies, noise 
peaks can help push objects that would not otherwise be 
detected over the detection threshold. (See also ^^4.31 ) 

Note that the above test explicitly avoids incomplete- 
ness due to source confusion. If we repeat the above test 
without avoiding known sources, we find that where com- 
pleteness is low, confusion actually increases the com- 
pleteness by a factor of a few, with faint sources hid- 
ing i n the skirts of brighter ones (see also iBerta et al.l 
l2006f) . However, the flux measurements for these ob- 
jects are naturally dominated by their neighbours; in this 
sense, it is arguable as to whether the synthetic object 
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Fig. 4. — Completeness for synthetic i?^/''-law sources added 
to the I'.'O FWHM K detection mosaic — At each point, we give 
the fraction of synthetic sources (characterized by their total K 
magnitude and effective or halflight radius, -Reff) that we are able 
to recover, using identical processes as for 'live' detection. Each 
synthetic source has been isolated from any other real or syn- 
thetic source by at least 50 pix to prevent confusion. While we 
are 100 % complete for point sources (i.e. i?;,ff = 0) at our nominal 
limit of K = 22, the completeness drops rather rapidly for larger 
radii/lower surface brightnesses. At fixed K and R^ft, both sur- 
face brightness and completeness are strong functions of the profile 
shape; we present these results for R^^'^ sources as loose lower lim- 
its on the true completeness. 
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Fig. 5. — Spatial Variation in the completeness of the MUSYC 
ECDFS catalogs — Completeness for synthetic K = 22.4 point 
sources randomly added to the K detection image. Completeness 
is slightly lower around the very bright star towards the eastern 
edge of the field, but is otherwise reasonably uniform. 



is actually being 'detected'. Where completeness is high 
{K < 20.5), confusion reduces completeness by a few 
percent, but again, the exact amount is sensitive to the 
position and flux agreement required to define a success- 
ful detection. From these tests, it seems that < 2 % of 
sources are affected by confusion due to chance align- 
ments with foreground/background galaxies (c/. gravita- 
tional associations). For comparison, based on the SEx- 
tractor segmentation map, if-detected objects cover 2.34 
% of the field. 

We have also done a similar test to investigate any vari- 
ations in completeness across the field. We placed 5000 
point sources with K — 22.4 — our 50 % completeness 
limit for point sources — across the field, each isolated by 
at least 26'.'7 (100 pix). The results are shown in Figure 
O Although it is perhaps slightly lower for the noisier 
east and northeast pointings, the completeness is indeed 
quite uniform across the full field. 

Finally, we can obtain empirical measures of both com- 
pleteness and reliability by comparing our catalog to the 
much deeper j^-select ed FIREWORKS catalog of the 
GOODS-CDFS region (iWuvts et al. l l200l The results of 
this exercise are shown in Figure [SI Here, the complete- 
ness is just the fraction of FIREWORKS sources which 
also appear in the MUSYC catalog; similarly, the reliabil- 
ity is the fraction of MUSYC sources which do not appear 
in the FIREWORKS catalog. For the 21.8 < K < 22.0 
bin, the MUSYC catalog is 87.5 % complete, and 97 % 
reliable. For K < 22, the overall completeness and relia- 
bility are 97 % and 99 %, respectively. 

Since the GOODS-ISAAC data are so much deeper, 
the high completeness at K 22 implies that K ~ 22, 
Re ^ 0'.'5 objects make up at most a small fraction of 
the FIREWORKS catalog. This might imply that our 
catalog is primarily fiux, rather than surface brightness, 
limited. It must also be remembered, however, that the 
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Fig. 6.— Completeness and reliability of the MUSYC ECDFS 
catal og in comparison to the FIREWORKS catalog (Wuv ts et al.l 
[2003) of the GOODS-CDFS data— We show the fraction of FIRE- 
WORKS sources that are detected by MUSYC (i.e., the complete- 
ness of the MUSYC catalog; solid black histograms), and the frac- 
tion of MUSYC sources that are confirmed by FIREWORKS (i.e., 
the reliability of the MUSYC catalog; dashed red histograms), in 
bins of total apparent K magnitude. For the 21.75 < K < 22.00 
bin, we are more than 85 % complete, with essentially all detections 
confirmed by the (much deeper) GOODS data. Cumulatively, to 
K < 22, the MUSYC catalog is 97 % complete (black curve), and 
better than 99 % reliable (red curve). 

main motivation for large area surveys like MUSYC is to 
find the rare objects that may be missed in smaller area 
surveys hke GOODS. 

4.3. Total Fluxes — Method 

We measure total fluxes in the I'.'O FWHM K band 
mosaic, using SExtractor's AUTO measurement, which 
uses a flexible elliptical aperture whose size ultimately 
depends on the distribution of light in 'detection' pix- 
els (i.e., an isophotal region). We do specify a 
minimum AUTO aperture size (using the parameter 
PHOT_AUTOAPERS) of 2'.'5, although in practice this 
limit is almost never reached for sources with K < 22. 
The 2'.' 5 limit has been chosen to be small enough to 
ensure high signal-to-noise for faint point sources, while 
still avoiding any significant aperture matching effects 
(see both §3.21 and §4.6p . We apply two corrections to 
the AUTO flux to obtain better estimates of galaxies' 
total fluxes; these are described below. We will then 
quantify the effect and importance of these corrections 
in the following section. 

Even for a point source, any aperture that is 
comparable in size to the PSF will miss a non- 
negligible amount of flux (e.g. iBertin fc Arnoutd 
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1996 



2002 



Labbe et al.l 120031 : iBrown et all l2007f) . iBrown et ail 
()2007| ) have shown that fraction of light missed by the 
AUTO aperture correlates strongly with total magni- 
tude; this is simply due to the fact that the AUTO 
aperture siz e corr elates strongly with total brightness. 
iLabbe et al.l ()2003f ) fin d that up to . 7 mag can be missed 
for some objects, and IBrown et al.l ()2007D suggest that 
the systematic effect at the faint end is ^ 0.2 mag. 

It is therefore both appropriate and important to ap- 
ply a correction for mis sed fiux laying out side the 'to- 
tal' aperture. Following ILabbe et al.l (|2003t ). we do this 
treating every object as if it were a point source: us- 
ing the empirical models of the PSF constructed as per 
H3.21 we determine the fraction of light that falls outside 
each aperture as a function of its size and ellipticity, and 



scale SExtractor's FLUX_AUTO measurement accord- 
ingly. Since no object can have a growth curve which is 
steeper than a point source, this is a minimal correction: 
it leads to a lower limit on the total flux. 

Further, we flnd that SExtractor's background esti- 
mation algorithms systematically overestimate the back- 
ground level, which also produces a bias towards lower 
fluxes. Because SExtractor does not allow the user to 
turn off background subtraction when doing photometry 
(c/. detection), we are forced to undo SExtractor's back- 
ground subtraction for the final catalog, using the out- 
put BACKGROUND values, and the area of the AUTO 
aperture. We have done this only for the total K fiuxes; 
since we have performed our own background subtrac- 
tion for the NIR images (as described in i j2.2.3[) . undo- 
ing SExtractor's background subtraction is equivalent to 
trusting our own determination. Note that, for the SED 
fluxes, we still rely on SExtractor's LOCAL background 
subtraction algorithm, with PHOTO.THICK set to 48. 



4.4. Total Fluxes 



Validation 



Our overarching concern here is the correspondence be- 
tween our measured fluxes and the true total fluxes of real 
sources. We have tested our total flux measurements by 
checking our ability to recover the known fluxes of large 
numbers synthetic sources, inserted into the I'.'O FWHM 
K science image as in i |4.2l The results of these tests 
are shown in Figure [71 In this Figure, we compare the 
performance of SExtractor's AUTO measurement before 
(upper panels) and after [lower panels) our corrections 
for missed flux and background over-subtraction are ap- 
plied. In each case, the contours show the systematic 
(left panels) and random [right panels) errors in the re- 
covered magnitude. The red lines show the approximate 
90 %, 50 %, and 10 % completeness limits for R^^^-lim 
sources, as derived in §4.21 

Further, in order to gauge the way these measurements 
are affected by noise, we have performed several vari- 
ations of this test. In each test we add the synthetic 
sources either to a noiseless image, or to the actual I'.'O 
FWHM K mosaic; we have trialled the four possible per- 
mutations of using the noiseless or real image for detec- 
tion or measurement. We briefly summarize the results 
of these tests below. 

The reader wishing to avoid such a technical discussion 
of SExtractor's photometry algorithms may wish to skip 
to i i4.5l after noting that, comparing the upper and lower 
panels of FigurejTl the effect of our two corrections to the 
AUTO measurement is to reduce the systematic underes- 
timate of total fluxes by > 0.10 mag. For point sources, 
the total flux is recovered to within 0.02 mag for K < 22. 

4.4.1. Missed Flux and Aperture Size Effects 

In order to determine the bias inherent in the AUTO 
algorithm, we have checked our ability to recover the 
fluxes of synthetic sources placed in a noiseless image, 
using this image for both detection and measurement. 
For point sources, the photometric bias inherent in the 
AUTO algorithm is < 0.05 mag for K < 20.5, but rises 
to 0.10 mag for K = 22. It is also a strong function of 
RcS- aX K = 21.5, the AUTO aperture misses 0.12 mag 
for i?cff = 0'.'4, and more than 0.25 mag for i?cff — I'-'O. 
Applying our 'point source' correction for missed flux 
reduces this bias to < 0.02 mag for all K < 22 point 
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Fig. 7. — Validating our total flux measurements — The systematic {lejt panels) and random (right panels) errors in the recovered fluxes 
of synthetic i?^/*-law sources introduced into the I'.'O FWHM K science image, based on SExtractor's AUTO aperture, before {upper 
panels) and after {lower panels) applying corrections for missed flux and background oversubtraction. The red lines in each panel show 
the approximate 90 % (solid), 50 % {long dashed), and 10 % (short dashed) completeness limits for i?^/* law sources, as in Figure|4] As in 
Figure |4] the results shown in this Figure are presented as upper limits on the systematic errors; both the random and systematic errors are 
significantly less assuming exponential profiles. In order to account for flux laying beyond the AUTO aperture, we correct the flux of each 
source as if it were a point source; this is thus a minimal correction. This correction reduces the systematic error in total fluxes by ^0.1 
mag for an Rl/-* profile, and from < 0.10 mag to < 0.01 mag for point sources. For the faintest sources, this correction also reduces the 
random error by as much as 0.05 mag. We also find that SExtractor's LOCAL background subtraction algorithm tends to over-estimate 
and over-subtract the background, leading to a systematic offset at the level of 0.03 mag. With both of these corrections, the systematic 
errors in total fiuxes for point sources are reduced to < 2 % 



sources; and, at K — 21.5, to 0.08 and 0.21 mag for 
ReB = 0'.'4 and I'.'O, respectively. 

The above numbers indicate the bias inherent in the 
AUTO algorithm, even for infinite signal-to-noise; con- 
sidering synthetic sources introduced into the real K 
science image, we find that noise exacerbates the prob- 
lem. For point sources, the mean offset between the un- 
corrected AUTO and total fluxes are < 0.05 mag for 
K < 20.0, 0.10 mag for K = 21.5 and 0.17 mag for 
K = 22.0. For K = 21.5, the systematic offset is 0.16 
mag for Res = 0'.'4 and 0.50 mag for Res = I'.'O. For 
K = 22, the average 'point source' correction for missed 
flux goes from 0.05 mag for true point sources up to 0.10 
mag for i?cff = I'-'O, and 0.15 mag for R^ff = 1'.'5. After 
applying our correction for missed flux, the photometric 



offset is reduced to < 0.03 mag for all point sources; at 
K = 21.5, the numbers for i?cff = 0'.'4 and ReS = I'.'O 
become 0.10 mag, and 0.35 mag, respectively. 

As an aside, we have also looked at how noise in the de- 
tection image affects the AUTO measurement, by using 
the real image (with synthetic sources added), for detec- 
tion, and using a noiseless image for measurement. The 
effect of noise in the detection image is to induce scatter 
in the isophotal area, and so the AUTO aperture size, at 
a flxed i?eff and K. Applying a correction for missed flux 
thus reduces the random scatter in the recovered fluxes 
of low surface brightness sources, by eliminating the first 
order effects due to aperture size; the random scatter in 
recovered fluxes is reduced by ~ 0.05 mag for all K <21 
sources. This can be seen in Figure [T] 
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Also, as in note that the numbers given above 

all apply to galaxies with an i?^/"* profile, and so should 
be treated as approximate upper limits on the random 
and systematic errors. We have performed the same test 
assuming exponential profiles: the systematic error in 
the recovered flux is less than 0.03 mag for all K < 22 
and i?cff < 0'.'6. 

4.4.2. Background Oversubtraction 

Even after correcting for missed flux, and even for point 
sources, SExtractor's photometry systematically under- 
estimates the total fluxes of synthetic sources. At least 
part of this lingering offset is a product of the LOCAL 
background subtraction algorithm. This algorithm uses a 
'rectangular annulus' with a user-specified thickness, sur- 
rounding the quasi-isophotal detection region. Any flux 
from the source lying beyond this 'aperture' (which may 
well be smaller than the AUTO aperture) will therefore 
bias the background estimate upwards, leading to over- 
subtraction, and so a systematic underestimate of the 
total flux. 

If we undo SExtractor's background subtraction^, then 
the photometric offset for point sources is reduced to 
< 0.02 mag for all K < 22. The size of this correction is 
only weakly dependent on source size and flux, varying 
from > -0.025 mag for {K, i?eff) = (19, 0'.'4) to -0.038 
for (K, R,s) = (22, 0'.'4). 

4.5. Multi-color SEDs 

In order to maximize signal-to-noise for the faintest 
objects, instead of measuring total fluxes in all bands, 
we construct multi-color SEDs based on smaller, 'color' 
apertures; we then use the K band total flux to normalize 
each SED. 

The 'color' photometry is measured from 1'.'5 FWHM 
PSF-matched images (see ij3.2p . again using the I'.'O 
FWHM K mosaic as the detection image. Specifically, 
we use SExtractor's MAGJSO, again enforcing a mini- 
mum aperture size of 2'.'5 diameter. This limit is reached 
by essentially all objects with K > 21.5, and essentially 
none with K < 20.5. Note that, even though the ISO 
aperture is defined from I'.'O FWHM K mosaic, (after 
SExtractor's internal filtering; see ij4.ip . all 'color' mea- 
surements are indeed made using matched apertures on 
1'.'5 FWHM PSF-matched images. 

In order to test our sensitivity to color gradients, we 
have verified that i?tot = {R — K)co\ + -K^tot, where Rtot 
comes from using the R band image in place of the K 
band image for detection and total flux measurement. 
Particularly for the brightest and biggest (and so, pre- 
sumably, the nearest) galaxies, the use of the ISO aper- 
ture is crucial in ensuring that this is indeed true. 



we empirically determine the photometric measurement 
uncertainties by placing large numbers of apertures on 
empty or blank regions in our measurement images. 
The principal advantage of this approach is that it 
correctly accounts for pixel-pixel correlations intro- 
duced in various stages of the data reduction process 
(including interpolation during astrometric correction 
and convolution during PSF matching). 

For the 'color' apertures, we have placed 2'.'5 — 8" inde- 
pendent (i.e., non-overlapping) apertures at 10^ 'empty' 
locations, based on the combined optical {BVR) and 
NIR [K) segmentation maps. With this information, 
we can build curves of <7{A) for each band, where a is 
the measurement uncertainty in an aperture with area A. 
Similarly, for the 'total' apertures, which are somewhat 
larger, we have placed 2'.'5-12" independent apertures at 
3500 'empty' locations on the I'.'O FWHM K detection 
mosaic, using only the NIR segmentation map to define 
'empty'. Note that since the 'empty aperture' photome- 
try is done using SExtractor, in the same manner as for 
our final photometry, the errors so derived also account 
for random uncertainties due to, for example, errors in 
background estimation, etc. 

There is one additional layer of complexity for the ISPI 
bands: in order to track the spatial variations in the 
'background' RMS, both within and between subfields, 
we use the RMS maps produced during mosaicking by 
xdimsum (see ij2.2.2p . While these maps are not ac- 
curate in an absolute sense, they do adequately map 
the shape of RMS variations across each subfield. We 
have therefore normalized these maps by the RMS flux 
in empty 2'.'5 apertures, and then combined them to con- 
struct a (re)normalized 'RMS map' for the full 30' x 30' 
fleld. Then, in practice, the photometric uncertainty for 
a given object is estimated by taking the median pixel 
value within the SExtractor segmentation region associ- 
ated with that object, corrected up from 2'.'5 to the ap- 
propriate aperture size using the (j{A) curves described 
above. 
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4.6. Photometric Errors 
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^ Again, note that SExtractor does not allow the user to turn 
off background subtraction for photometry. In practice, we have 
undone SExtractor's background subtraction using the output 
BACKGROUND value, multiplied by the area of the AUTO aper- 
ture. The AUTO aperture area is given by KRON_RADIUS^ X tt X 
AJMAGE X BJMAGE. Note, too, that we apply this correction 
before the missed flux correction discussed in i|4.4. II 



Fig. 8. — Validating our error estimates — For the JK imaging, 
we use the RMS maps output by xdimsum to account for spatial 
variations in the background noise level. Although these maps are 
not accurate in an absolute sense, they do adequately map the 
relative variations in the noise; accordingly, we have normalized 
these RMS maps using the scatter in 2'.' 5 diameter apertures placed 
on empty regions of the science image. The points in this Figure 
show the integrated flux in each of 3000 of these apertures, F, as 
a function of the estimated error, AF, derived using these RMS 
maps. The red error bars show the RMS in F, binned by AF. The 
observed RMS in empty apertures agrees extremely well with the 
estimated errors based on the normalized xdimsum RMS map. 
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In Figure m we validate these error estimates by show- 
ing the 'empty aperture' fluxes, F, as measured in 2'.'5 
diameter apertures, as a function of the photometric er- 
ror, A.F, estimated as above. The line with error bars 
shows the mean and variance of the 'empty aperture' 
fluxes in bins of AF; in other words, the error bars show 
the actual error, plotted as a function of the estimated 
error. The agreement between the photometric errors es- 
timated using the RMS map, and the variance in 'empty 
aperture' fluxes is excellent. This is more than just a 
consistency check: while the RMS maps have been nor- 
malised to match the variance in empty aperture fluxes 
on average, the fact that the observed scatter scales so 
well with the predicted error demonstrates that the RMS 
map does a good job of reproducing the spatial variations 
in the noise. 

For a Gaussian profile {i.e., a point source), and in 
the case of uncorrelated noise, an aperture with a di- 
ameter 1.35 times th e FWHM gives the optimal S:N 
(jGawiser et al.l [2006al) . Based on the 'empty aperture' 
analysis described in ij4.61 the 2'.'5 aperture size is slightly 
larger than optimal for a point source in the J (I'.'S 
FWHM) image. For the I'.'O FWHM K detection image, 
the optimal aperture diameter for a point source is 1'.'33; 
the S:N in a 2'.'5 diameter aperture is 25 % lower. Us- 
ing slightly larger apertures presumably increases S:N for 
slightly extended sources, as well as reducing sensitivity 
to systematic effects due to various classes of aperture 
effects (e.g., imperfect astrometric and PSF matching, 
etc.). 

Within a 2'.'5 diameter aperture, the formal ba limits in 
the K band are 22.25 mag at an effective weight of 0.75, 
and 22.50 mag at an effective weight of 1.0. Averaged 
across the image, the 5(7 limit is 22.42 mag; the limits 
for all bands are given in Table [TJ For a point source, 
these limits can be translated to total fluxes by simply 
subtracting 0.45 mag. 

5. ADDITIONAL CHECKS ON THE MUSYC 
CALIBRATION 

5.1. Checks on the Astrometric Calibration 

In order to test the astrometric calibration of the 
MUSYC ECDFS imaging, we have compared the cat- 
aloged position of sources from the if -selected catalog 
with those from version 3.3 of the Y ale/San Juan South - 
ern Proper Motion (SPM) catalog (|Girard et al.ll2004D . 
This catalog is based on observations made using the 
51 cm double astrograph of Cesco Observatory in El 
Leoncito, Argentina. For V < 17, the positional accu- 
racy of the catalog is 0'.'04— 0'.'06. 

In Figure[51 we show an astrometric comparison for 113 
objects common to the SPM and MUSYC catalogs; these 
objects are plotted as black circles. For this comparison 
we have selected objects with 14 < U < 18 and proper 
motions of less than 20 mas / year. All these objects 
have 14 < if < 18; the median has K — IQ mag. 

The systematic offset between SPM- and MUSYC- 
measured positions, averaged across the entire field, is 
0'.'079 in Right Ascension and 0'.'222 in declination; that 
is, a mean offset of 0'.'235 (0.88 pix), 20° East of North. 
For these sources, the random error in the MUSYC po- 
sitions is 0'.'30 and 0'.'27 in x and y, respectively. 

We have performed the same comparison for the 
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Fig. 9.— Validating the MUSYC ECDFS astrometric 
calibration — In the upper part of this Figure, we show the field- 
averaged astrometric comparison between the MUSYC ECDFS 
catalog, and the Yale/San Juan Southern Proper Motion (SPM) 
catalog v3.3 (circles), as well 2MASS (crosses). In the flank- 
ing panels, the solid (dashed) histograms show the distribution of 
RA/dec offsets with respect to the SPM (2MASS) catalogs; we also 
give the mean and RMS offset between MUSYC and SPM catalog 
positions. In the lower part of this Figure, we show astrometric off- 
sets as a function of position; in these panels, the solid grey lines 
show the median-filtered relation derived from the SPM points. In 
comparison to the SPM catalog, the MUSYC astrometry is offset 
by 0'.'23 (0.87 pix); there is also evidence of an astrometric shear 
of < O'.'S (1.1 pix) in the RA direction across the full field. 



2MASS sources that were used in the photometric cal- 
ibration of the K images; these objects are shown in 
Figure 2 as the grey crosses. The median K magni- 
tude of these objects is 14.75 mag, considerably brighter 
than the SPM sources used above. In comparison to the 
2MASS catalog, which has astrometric accuracy of < O'.'l 
for K < 14, we find a slightly larger random offset: (0'.'22, 
0'.'39) in (RA, dec). For these sources, the random error 
in (RA, dec) is (0'.'22, C'lg). 

In the lower part of FigurelHl we plot the positional off- 
sets as a function of position accross the field. In these 
panels, the solid grey line shows the median- filtered re- 
lation between SPM- and MUSYC-measured positions. 
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There appears to be a slight astrometric shear in the RA 
direction at the < O'.'S level from the East to the West 
edge of the K mosaic. Otherwise, however, the offsets are 
consistent with the direct shift of 0'.'235 derived above. 

5.2. Checks on the Photometric Calibration 
5.2.1. Comparison with FIRE WORKS 

In order to test our photometric calibration, we have 
compared ou r catalog to the FIREWORKS catalog 
(|Wuvts et all HOOS) of the GOODS-CDFS region (the 
central ~ 150 □" of our field), which includes HST-ACS 
optical imaging, and significantly deeper NIR imaging 
taken using ISAAC on the VLT. Since the FIREWORKS 
catalog uses different filters, we are forced to use stellar 
colors to make this comparison. The results of this com- 
parison are shown in Figure fTUl Each panel in this Figure 
shows the color-color diagram for stars in terms of their 
FIREWORKS {Veoew - hrsw) color, and a MUSYC- 
minus-FIREWORKS 'color'. In each panel, the circles 
with error bars show the observations; these error bars 
apply only to errors in the MUSYC photometry. 

We have used spectra for lumin osity class V stars from 
the BPGS stellar spectral atlas (jGunn fc StrvkeHll983D 
to generate predictions for where the stellar sequence 
should lie in these diagrams. These predictions are the 
solid red lines in each panel; the small blue stars show the 
predicted photometry for individual BPGS stars. Note 
that, for the purposes of this comparison, we have con- 
verted to the Vega magnitude system, so that the stellar 
sequence necessarily passes through the point (0, 0). 

We calculate the photometric offset in each band as 
the S:N-weighted mean difference between the observed 
stellar photometry and the predicted stellar sequence. 
These values are given in each panel; the dashed red 
line is just the predicted stellar sequence offset by this 
amount. Our results do not change if we use the lPicklesI 
([1998) steUar atlas. 

Particularly for the NIR data, the absolute calibration 
of the MUSYC and FIREWORKS data agree very well: 
typically to better than 0.03 mag. In terms of the relative 
calibration across different bands, we see a discrepancy 
between the / and z' band calibrations of A(/— z') — 0.05 
mag, as well as a discrepancy between the U^s and B 
bands at the level of A{U3s — B) = —0.07 mag. 

5.2.2. Comparison with COMB 0-1 7 

Although the COMBO-17 broadband UsaBRVI imag- 
ing is a subset of the raw data used to produce the 
MUSYC imaging, the data reduction and analysis strate- 
gies used by each team are very different. For example, 
rather than a single measurement from a coadded im- 
age, the COMBO-17 flux measurements are based on 
the coadding of many distinct measurements from the 
individual exposures, and SED or 'color' measurements 
are made using adaptive, weighted 'apertures', rather 
than traditional (top-hat) apertures. Direct, object-by- 
object comparison between the two catalogs thus offers 
the chance to test both the photometric calibration, and 
the methods used for obtaining photometry. 

The results of this comparison are shown in middle 
panels of Figure [Tl] these panels show the difference 
in the MUSYC and COMBO-17 cataloged UssBVRI 
fluxes, plotted as a function of total R magnitude in 
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Fig. 10. — Photometric comparison between GOODS and 
MUSYC in the C DFS. based on the FIREWORKS catalog 
llWuvts et ahll2008l) of the GOODS data —In effect, these pan- 
els show the empirical color-transforms for stars between GOODS 
(ACS/ISAAC) and MUSYC (WFI/Mosaic-II/Son/ISPI) filters, 
plotted as a function of ACS color. Note that for this compar- 
ison, we have adopted the Vega magnitude system, so that the 
stellar sequence must necessarily pass through the point (0, 0). In 
each panel, the open circles with error bars represent the data; 
the errors shown here pertain only to the MUSYC photometry. 
The closed blue stars sh ow predicted stellar photometry based on 
the BPGS IIGunn &: Strvker.l983i^ stellar spectral atlas, convolved 
with the known filter curves; the solid red line shows a prediction 
for the stellar sequence in these color-color diagrams, obtained by 
median filtering the BPGS points. For each filter, we derive a pho- 
tometric offset by taking the mean difference, weighted by S:N in 
the MUSYC catalog, between the observed and predicted location 
of the stellar sequence in the y direction. These values are given in 
each panel (the 'A' having the sense of MUSYC— minus-GOODS); 
the dotted lines show the predicted stellar sequence offset by this 
amount. Particularly for the reddest bands, this comparison vali- 
dates the MUSYC photometric calibration at the few percent level. 
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Fig. 11. — Photometric comparison between COMBO-17 and 
MUSYC in the ECDFS — Each panel shows the difference be- 
tween the original COMBO-17 and MUSYC photometry; the red 
stars represent the observed photometry for stars, selected using 
COMBO-17's SED classification. In the central panels, which show 
direct object— by-object comparisons for those filters common to 
both MUSYC and COMBO-17, the black points show the same 
information for galaxies, plotted as a function of apparent R mag- 
nitude in the COMBO-17 catalog. The top and bottom panels 
compare the MUSYC U and z' bands to the COMBO-17 Uss and 
I bands, as a function of the MUSYC {B — V) color (in the Vega 
system), and based on synthetic photometry for main sequence 
stars. The median MUSYC-minus-COMBO-17 photometric off- 
set for each band is given in each panel; the dotted lines in each 
panel show the expected location of the stellar locus offset by this 
amount. Note that while the COMBO-17 raw data is a subset of 
the MUSYC raw data, the data reduction and analysis pipelines 
are completely independent. There are significant differences be- 
tween the COMBO-17 and MUSYC photometry, due at least in 
part to photometric cahb ration errors in the original COMBO-17 
catalog l lWolf et al][2003) . Even after recal ibrating the COMBO- 
17 photometry following fWolf et al.l'200i^, however, significant 
differences remain: for UgsBV RI , the offsets are —0.014, —0.141, 
-0.109, -0.112, and -0.124 mag, respectively. 



the COMBO-17 catalog, Ren- The comparison is be- 
tween total fluxes: i.e., /ci7 — Rtot.cn + {I ^ R)ci7] 
Imvs = KtotMUS + {I ~ K)mus- We have also trans- 
formed our data to the Vega magnitude system. For 
the purposes of this comparison, we distinguish between 
stars {red stars) and galaxies {black points), on the basis 
of the COMBO-17 SED classification; the results do not 
change significantly using Bz'K selected stars or GEMS 
point sources. We have used those stars with i? < 21 to 
identify differences in the two surveys' calibrations; these 
offsets are given in each panel, and shown as the dotted 
black lines. 

There are significant differences between the MUSYC 
and original COMBO-17 calibrations. These are due to 
calibration errors in the COMBO-17 catalog (Wol f et "all 
[2008). The original COMBO-17 calibration was based 
on spectrophotometric observations of two stars, each of 
which suggested different calibrations; in the end, the 
wrong star was chosen. Partially motiv ated by the 
comparison in Figure m IWolf et al l (2008) have since 
revised the basic calibration of the COMBO-17 ECDFS 
data using the other spectrophometric star, shifting the 
UssBVRI calibration by -0.143, +0M0, -f 0.003, -0.054, 
and -0.123 mag, respectively. 

We note that these rather large calibration errors do 
not have a huge effect on the COMBO-17 redshift deter- 
minations (jWolf et al.|[2008l . Paper II). This is because 
the medium bands, which are key to measuring break 
strengths and so choosing the redshift, are calibrated 
with respect to the nearest broad band. However, we 
show in Paper II that the effect on derived quantities 
like restframe colors and stellar masses is large. 

After recalibration using the other spectrophotomet- 
ric standard, the MUSYC and COMBO-17 stellar colors 
agree at the level of a few hundredths of a magnitude for 
BVRI; for U^s a discrepancy remains at the 0.1 mag 
level. Moreover, a discrepancy in the overall calibra- 
tion remains, such that stars are 0.1 mag brighter in the 
MUSYC catalog. Our correction for missed flux accounts 
for 0.03 mag of this offset; the source of the remaining 
0.07 mag offset has not been identified. 

Secondly, notice that there are apparently different off- 
sets for galaxies and stars: even after matching the two 
surveys' calibrations for stars using Figure llll galax- 
ies are still fainter and bluer in the COMBO-17 cata- 
log than they are in ours. Quantitatively, the UssBVRI 
galaxy-minus-star offsets are 0.102, 0.020, 0.010, 0.067, 
and 0.088 mag, respectively. Further, excepting the C/38 
band, the random scatter between the COMBO-17 and 
MUSYC galaxy photometry is 2 — 3 times greater than 
that for stars. It is difficult to say what might produce 
this effect, but the effect persists even when we use our 
R band image for detection and measurement; that is, 
this is not a product of our measuring total ffuxes in 
K rather than R. We do not believe that the combi- 
nation of COMBO-17's smaller effective apertures and 
galaxy color gradients can fully account for these effects. 
For R > 21, the effective diameter of the ISO aper- 
ture is almost always smaller than 2'.'5; for these objects 
the MUSYC photometry effectively uses fixed apertures. 

^" Note that these calibration issues affect only the ECDFS, and 
not the other three COMB O-17 fi elds, where multiple calibration 



stars give consistent results (|Wolf et al.ll2008l ). 
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TABLE 3 

Checks On the Photometric Calibration 



Ba,iid 


Photometric Offset with respect to 




FIREWORKS 


COMBO-17 


Steiiar SEDs 


CI ~1 








U 


+0.013 


+0.02 


-0.004 


Uas 


-0.020 


-0.15 


-0.051 


B 


+0.050 


-0.09 


-0.017 


V 


+0.038 


-0.09 


-0.006 


R 


+0.016 


-0.15 


+0.017 


I 


+0.055 


-0.23 


+0.023 


z' 


-0.004 


-0.27 


-0.011 


J 


+0.015 




+0.032 


H 


-0.012 




-0.032 


K 


-0.017 







Note. — This Table summarizes the results of i|5.2l For each 
band (Col. 1), we give: (Col 2) : the photometric offset between 
MUSYC and the FIREWORKS l lWuvts et al.l2008l) catalogs of the 
GOODS ACS and ISAAC imaging data; (Co L 3): the photom etric 
offset between the MUSYC and COMBO-17 l lWolf et al.ll2004l) op- 
tical imaging data; (Col. 4): the residuals from fitting stellar SEDs 
from the MUSYC catalog using main sequence stellar spectra from 
the BPGS atlas. 

While the agreement between star and galaxy colors is 
noticeably better for i? < 21 using fixed 2'.'5 apertures to 
construct SEDs, it does not have a significant effect for 
i? > 21, where the problem is greatest. 

While we cannot directly compare our U and z' pho- 
tometry to COMBO-17, it is still possible to use stel- 
lar colors to check these bands, as we have done for the 
FIREWORKS catalog. This is shown in the top and bot- 
tom panels of Figure [TT] For the z' band, this analysis 
suggests a possible discrepancy between the MUSYC / 
and z' band calibrations of A(/ — 2:')mus — 0.03 mag. 
For the U band, however, it suggests a discrepancy of 
A([/ — J738)mus ~ 0.15 mag. While we have been unable 
to identify the cause of this offset, we note both that the 
shape of the observed and predicted stellar sequences do 
not obviously agree as well for the U band as for the 
z' , and also that the results of both ii5.2.1l and i|5.2.3l 
do not support the notion of an offset of this size. We 
do not believe that this indicates an inconsistency in the 
calibrations of the U and bands. 

5.2.3. Refining the Photometric Cross-Calibration using 
Stellar SEDs 

In the construction of SEDs covering a broad wave- 
length range, the relative or cross-calibration across all 
bands is at least as important as the absolute calibration 
of each individual band. As a trivial example, if the 
zeropoints of two adjacent bands are out by a few 
percent, but in opposite senses, this can easily introduce 
systematic offsets in color on the order of 0.1 mag; 
the worry is then that these apparent 'breaks' might 
seriously affect photometric redshift determinations. 
This is a particular concern in the case of the MUSYC 
ECDFS dataset, which incorporates data from four 
different instruments, each reduced and calibrated using 
quite different strategies. 

We have therefore taken steps to improve the pho- 
tometric cross-calibration of the MUSYC ECDFS data. 
The essential idea here is to take a set of objects whose 
SEDs are known a priori (at least in a statistical sense) 



and to ensure agreement between the observed and ex- 
pected SEDs. Stars are, in fact, ideal for this purpose, 
since they form a narrow 'stellar sequence' when plotted 
in color-color space: at least in theory, and modulo the 
effects of, e.g., metalicity, a star's (c/. a galaxy's) full 
SED can be predicted on the basis of a single color. 

Our method is as follows. We begin with a set of 
more than 1000 objects with unambiguous 'Star' classi- 
fications in the COMBO-17 catalog, of which nearly 600 
have photometric S:N > 10 in K, and are unsaturated in 
all MUSYC bands. Again, our results do not change if 
we use Bz'K selected stars or GEMS point sources. Us- 
ing EAZY (see i j7.2l for a description), we fit the objects' 
photometry with luminosity class V stellar spectra from 
the BPGS stellar spectral atlas as a template set, and 
the redshift fixed to zero. Note that, by default, EAZY 
includes a 0.05 mag systematic error on each SED point, 
added in quadrature with the measurement uncertainty. 

Using the output to discard objects whose SEDs 
are not consistent with being a main sequence star, we 
can then interpret the median residual between the ob- 
served and best-fit photometry as being the product of 
calibration errors, and so refine the photometric calibra- 
tion of each band to ensure consistency across all bands. 
Specifically, given the photometric errors, we use min- 
imization to determine the zeropoint revision. 

The zeropoint revisions derived in this way are small; 
< 0.05 mag in all cases. The exact revisions are given in 
Table [21 Across the WFI data, there appears to be an 
offset that is roughly monotonic between the C/38 and R 
bands, where the offset in U^^ — I is —0.054 mag; cf. 0.055 
mag from the comparison to the FIREWORKS catalog. 
Similarly, there is an apparent inconsistency between the 
/ and z' calibrations, such that the offset in {I—z') is 0.03 
mag; cf. 0.05 mag from the comparison to FIREWORKS. 

The crux of this method is that whatever zeropoint 
discrepancies exist do not affect the choice of the best 
fit template in a systematic way. For example, a large 
offset in the U bands or a wavelength-dependent offset 
might lead to stars being fit with systematically bluer 
or redder template spectra, so biasing the derived pho- 
tometric offsets. In this sense, it is reassuring that the 
derived offsets are small, and comparable to the quoted 
uncertainties on the photometric calibration. Further, 
we note that we get very similar results if we increase 
the systematic uncertainty used by EAZY to 0.10 mag. 

Given the agreement between the results of the exter- 
nal comparison to FIREWORKS and those from the in- 
ternal consistency check on stellar colors, we have chosen 
to adopt the zeropoint revisions suggested by this stellar 
colors exercise. With these revisions, we believe that our 
photometric calibration is accurate, in both an absolute 
and a relative sense, to the level of a few hundredths of 
a magnitude. 

6. NUMBER COUNTS 

As a very basic comparison between our catalog and 
other /ir-selected catalogs. Figure [TH shows the number 
of detected galaxies as a function of total apparent K 
magnitude. Note that all the catalogs shown apply a 
similar correction for flux missed by SExtractor's AUTO 
measurement. The left panel of this figure shows the 
raw number counts; the right shows the number counts 
normalized by area. In both panels, it can be seen that 
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Fig. 12. — K band apparent magnitude number counts, comparing the MUSYC ECDFS catalog to other /sT-selected catalogs — Left 
Panel: the raw num bers of d etected sources, in bins of K magnitude for the ECDFS {heavy solid histogram), in comparison to: the other 
MUSYC wide fields l IBlanc et al. 200^, light solid h istog rams)\ the MU SYC deep fields (Quadri et al. 2007, light dashed hist ograms); the 
FIREWORKS catalog of the GOODS-CDFS data l|Wuvts et al.ll2008l . heavy dashed histogram); and the two FIRES fields ( Labbc et alj 
[2003; Forster-Schreiber ct al. 2006, red dashed histograms). Right panel: The normalized number counts for the same collection of datasets; 
the MUSYC ECDFS data are highlighted {heavy red points). At a fixed K magnitude, while the GOODS region of the ECDFS has 
approximately 80 % as many sources as are found in the ECDFS as a whole, in comparison to the other MUSYC wide fields, the ECDFS 
is underdense at the level of ^ 5 %. 



our number counts drop off for if > 22; our catalog is 
nearly, but not totally, complete for K = 22. 

The overall agreement between these different catalogs 
is very good. Assuming that the calibration of all cata- 
logs is solid, and looking at the left panel of Figure [T^l it 
can be seen that the ECDFS is slightly underdense — at 
the level of 4 — 6 % for 17.5 < K < 21.5. — in compar- 
ison to the two other MUSYC wide NIR selected cata- 
logs (jBlanc et al.ll2008f) . Conversely, the ECDFS number 
counts can be matched to the other two wide catalogs 
by adjusting the ECDFS K photometric calibration by 
-0.06 or -0.09 mag. 

In comparison to the number counts from the FIRE- 
WORKS catalog of the GOODS CDFS region, the 
GOODS region contains approximately 18 % fewer 
sources per unit area than the ECDFS as a whole. Even 
after matching the MUSYC ECDFS K calibration to 
the FIREWORKS catalog (see gSH), the GOODS re- 
gion remains underdense by 16 % in comparison to the 
ECDFS. 

7. PHOTOMETRIC REDSHIFTS 

7.1. Star/ Galaxy Separation 

We separate stars and galaxies from within the 
MUSYC ECDFS catalog on the basis of their Bz'K col- 
ors. The Bz'K diagram is known as a means of selectin g 
moderate redshift (z > 1.4) galaxies (|Daddi et al.ll2004 ). 
but can also be used as a efficient me ans of distin- 
guishing s tars fro m galaxies (see, e.g., 'Grazia n et al.l 
[2006; Blanc et all [20081 . In Figure [Tl we evaluate the 
performance of this criterion in comparison to the stellar 



SED classification from COMBO-17 (|Wolf et all 12004V 
as we ll as to a cata log of point sources from GEMS 
(Hau ssler et al.ll2007f) . 

Both panels of Figure [T3| show the Bz'K diagram for 
the MUSYC ECDFS catalog {Uack points)] the Bz'K 
stellar selection line: 

(z' - if) < 0.3 (S - z') - 0.5 , (1) 

is shown as the dashed line. In total, from the main 
K < 22 sample, 755 sources are selected as stars on the 
basis of their Bz'K colors. The left-hand panel of Figure 
[13| shows where Bz'K star selection agrees with other in- 
dicators; the right-hand panel shows where there is dis- 
agreement. For instance, on the left, the star-shaped 
symbols show objects that are classified as 'stars' by 
COMBO-17; on the right, they represent those Bz'K- 
selected 'stars' which are not classified as such in the 
COMBO-17 catalog. Similarly, the circles refer to point 
sources in the GEMS catalog. In both panels, objects 
that have been spectrally identified as stars are high- 
lighted in red. In either panel, the stellar sequence is 
immediately obvious and, for a given (B — z') color, can 
be seen to be separated from the galaxy population by 
at least a few tenths of a magnitude in {z' — K). 

Looking at the left panel, there is near complete over- 
lap between COMBO-17's star classification and Bz'K 
selection: only a very few COMBO-17 'stars' lie above 
the Bz'K selection line. There are a few dozen GEMS 
point sources found above the Bz'K selection line. In 
the MUSYC and GEMS optical images, some are clearly 
non-circular, and only a few show diffraction spikes; 
these appear to be compact, un- or barely-resolved 
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Fig. 13. — Stellar identification using Bz'K colors — In each panel, we show the Bz' K diagram for sources in the MUSYC ECDFS 
catalog (black points), and compare our Bz' K star selection {dashed line) to other compl imentary stellar cla ssifications: viz. SED-classified 
'stars' from the COMBO-17 survey ( Wolf et al., .2004. open stars), GEMS point sources IIHaussler et al.ll20071 . open circles), and spectrally 
classified stars {open squares/red symbols). The left panel shows the agreement between Bz' K selection and these other indicators; in the 
right panel we show where Bz'A'-selection disagrees with other indicators. So, for example, circles in the left panel show all GEMS point 
sources, whereas in the right panel they show those iJ^'/C -selected 'stars' that are not GEMS point sources. In the either panel, the stellar 
sequence in Bz' K color space can be seen to be isolated by > 0.1 mag in {z' — K) from deep field galaxies. This includes QSOs, which can 
be seen in the left panel as GEMS point sources scattered throughout the galaxy population. Although there are a handful of spectrally 
classified stars lying well outside the Bz' K stellar selection region (open squares in the left panel), these objects are neither COMBO-17 
'stars' nor GEMS point sources (stars and circles in the right panel); i.e. the spectral classification is wrong. Of the Bz'ii'— selected stars 
which are not GEMS point sources (circles in the right panel), roughly half are faint stars superposed over a diffuse background galaxy, 
and roughly half are faint galaxies whose photometry is significantly affected by a bright, nearby star. 



galaxies. Note, too, that this region of the Bz'K dia- 
gram is sparsely populated by X-ray s ources {i.e. QSOs; 
iDaddi et alll20'o4 iGrazian et al.ll2006D . 

There are also a handful of objects that are spectro- 
scopically identified as stars, which also fall above the 
Bz'K star selection line. With one exception, however, 
these objects are not GEMS point sources (squares in the 
left panel; circles in the right); neither are they classified 
as stars by COMBO-17 (squares in the left panel; stars in 
the right). These are, therefore, probably erroneous spec- 
tral classifications. There are no spectroscopic galaxies 
that lie in the stellar region of the Bz' K diagram. 

Turning now to the right panel, there are 66 Bz' K- 
selected 'stars' which do not appear in the GEMS point 
source catalog. A handful of these simply did not receive 
GEMS coverage. Of the rest, visual inspection shows 
these sources to be, in roughly equal proportions, faint 
stars superposed over a faint, background disk galaxy, or 
faint galaxies whose photometry is significantly affected 
by a nearby bright star. There are also 76 Bz' K-seXecteA 
'stars' which are not classified as such in the COMBO- 
17 catalog. In (J — K)-K color-magnitude space, these 
objects almost all have (J — K) < and K < 21; this 
would suggest that these are faint stars misclassified by 
COMBO-17. 



7.2. Photometric Redshifts — Method 



The basic idea behind photometric redshift estimation 
is to use the observed SED to determine the proba- 
bility of an object's having a particular spectral type, 
t (drawn or constructed from a library of template 
spectra), and being at a particular redshift, z: i.e. 
p{z, i|SED). We have derived photometric redshifts for 
every object in the catalog using a new photometric red- 
shift code called EAZY (Easy and Accurate 2:photS from 
Yale; for a more det ailed and complete discussion, see 
iBrammer et al.l [2008[) . EAZY combines many features 
of other commonly used photometric redshift codes like 
a Bayesian luminosity prio r (e.g. BPZ; ,Benftez 2000[) 
and template combination (|Rudnick et al.l 1200 iL 120031 ) 
with a si mple user interface based on the popular hy- 
perz code (jBolzonella. Miralles k Pell6ll2000f) . Novel fea- 
tures include the inclusion of a 'template error func- 
tion'; a restframe wavelength dependent systematic er- 
ror, which down- weights those parts of the spectrum like 
the restframe UV, where galaxies show significant scat- 
ter in color-color space. Moreover, the user is offered 
full control over whether and how these features are em- 
ployed. 

Another key difference is that objects are assigned 
redshifts by taking a probability weighted integral 
over the full redshift grid {i.e. marginalizing over the 
posterior redshift probability distribution), rather than, 
for example, choosing the single most likely redshift. 
(Although again the user is given the choice of which 
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TABLE 4 

Summary of the Contents of the Photometric Redshift Catalog 



Column No. Column Title Description 



1 


id 


Object identifier, beginning from 1, as in the photometric catalog 


2 


z_spcc 


Spectroscopic redshift determination, where available, as given in the photometric catalog 


3, 4 


z_a, chi_a 


Maximum likelihood redshift, allowing non-negative combinations of all six of the default 
EAZY templates, and the value associated with each fit 


5, 6 


z_p, chi_p 


As above, but with the inclusion of a iC luminosity prior 


7,8 


z_ml, z_m2 


Probability-weighted mean redshift, without and with the inclusion o{ a K luminosity prior, 
respectively; we recommend the use of the z_m2 redshift estimator. 


9—14 


168, u68, etc. 


Lower and upper limits on the redshift at 68, 95, and 99 % confidence, as computed from 
the same posterior probability distribution used to calculate z_m2 


15 


odds 


The fraction of the total integrated probability within ±0.2 of the z_m2 value 


16 


qz 


The figure of merit proposed bvlBrammer et al.l dSOOSA, calculated for the z_m2 value 


17 


nfilt 


The number of photometric points used to calculate all of the above 




0.0 0.5 1.0 1.5 2.02.53.0 0.0 0.5 1.0 1.5 2.02.53.0 

COMBO- 17 photometric redshift GOODS-MUSIC photometric redshift 



Fig. 14. — Validating the MUSYC ECDFS photometric redshifts — Each panel shows an object-by— object comparison between the 
MUSYC photometric redshift, and that from COMBO-17 (Wolf ct al. 2004, left panel), and the GOODS-MUSIC catalog of the GOODS- 
CDFS data IGrazian et al..,2006 . right panel). In order to discriminate between the two ^phot^ where there is disagreement, where a robust 
spectroscopic redshift determination is available (see Appendix |A]| i the red lines connect each point in (zphoti ^phot) point to the point 
(zspcc, 2spGc); vertical lines thus indicate catastrophic errors in the MUSYC Zphots, where horizontal lines show catastrophic failures in the 
COMBO-17/GOODS-MUSIC .SphotS. COMBO-17 suffers from a few different classes of systematic effects, owing principally to the lack 
of NIR data. Note, however, that very few spectroscopic redshifts are available for these objects — these effects would not be noticeable in 
a Zgpec — Zpijot diagram. In this panel. X-ray selected sources are marked with a cross; within this R < 24 sample. X-ray-selected sources 
are roughly three times as likely to have |Az|/(l + ^spoc) > 0.1 (see also Appendix Rt. In comparison to GOODS-MUSIC, the MUSYC 
Zphot^ have a slightly greater number of catastrophic outliers, such that the MUSYC Zphot too low; again, many of these objects are 

X-ray sources. In this panel, objects with poorly constrained photometric redshifts (i.e. Qz > 3) are marked with a circle; these objects 
are roughly twice as likely to have |(zmuS ~ -^003)1/(1 + -^GDs) > 0-2. The overall agreement between the two redshift determinations is 
really very good, especially moving towards the 'redshift desert' at Zphot ^ 1-5. 
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estimator to use.) EAZY also outputs 68/95/99 % 
confidence intervals, as derived from the typically 
asymmetric p{z). EAZY thus outputs meaningful and 
reliable photometric redshift errors, including the effects 
of 'template mismatch'; i.e. degeneracies between the 
redshift solution and the spectral type. By Monte 
Carlo'ing our catalog (i.e. reanalyzing many Monte 
Carlo realizations of our photometry, perturbed accord- 
ing to the photometric errors), we have verified that the 
EAZY p{z) does in fact provide a good description of 
the redshift uncertainties due to photometric errors. 

We have adopted EAZY's default parameter set for 
our redshift calculations.^^ That is, we use a library of 
six template spectra, allowing non-negative linear com- 
binations between these basis templates, and including 
an apparent K magnitude prior, p(z\K), and using the 
default EAZY template error function. 

Both the base t emplate set an d the K prior have 
been derived by iBrammer et al.l (|2008l ) using syn- 
thetic pho tometry from the semi-analytic model of 
iDe Lucia fc Blaizot| (I2007D. which is in turn based on the 
Millenium Simulation ( Springel et al.|[2005f ). The moti- 
vation for this approach is to approximately account for 
the full diversity in < 2; < 4 galaxies' SEDs due to dif- 
ferences in their individual star formation and assembly 
histories. The K p r ior is constructed directly from the 
iDe Lucia fc BlaizotI (|2Q03) simulation. 

In order t o derive the base template set, 

IBrammer et al.l (120081 ) have applied th e non-negative ma- 
trix fa ctorization (NMF) algorithm of lBlanton fc Roweid 
(|2007f ). to this synthetic catalog. In essence, this algo- 
rithm takes a large template library and distills from it 
a reduced set of basis templates that best describe the 
full range of 'o b served ' photometry. For this purpose, 
IBrammer et al.l (120 08") have used the template library 
used by iGrazianet al. (2000) to generate photometric 
redshifts for the GOODS-MUSIC catalog. This library 
consists of ~ 3000 Pegase synthetic spectra with a 
variety of dust obscuration, star formation histories, and 
ages. In additional to the f ive base t e mplat es output by 
the NMF algorithm, Bram mer et al.l (j2008f ) also include 
one young, dusty template {t = 50 Myr; Ay — 2.75), 
to compensate for the l ack of dusty galaxies in the 
IDe Lucia fc BlaizotI (l2007f) similuation. 

iGrazian et all (|2006t r using their full template library, 
achieved a photometric redshift accuracy of Cz = 0.045 
for their GOODS-MUSIC catalog of the GOODS ACS- 
ISAAC-IRAC data. For the same data, and using the 
default setup described above, the EAZY photometric 
accuracy is ct^ = 0.036. This represents the current state 
of the art for photometric redshift calculations based on 
broadband photometry. 

Table |4] gives a summary of the information contained 
within the photometric redshift catalog. Note that when 
computing photometric redshifts, we only use photome- 
try with an effective weight of 0.6 or greater. In addition 
to the basic EAZY output, we have included two addi- 
tional pieces of information. The first is simply a binary 
flag indicating whether or not each object is classified as 

In Paper II, we present a number of variations on the pho- 
tometric redshift computation described here; in relation to Paper 
II, the redshifts described here correspond to the 'default analysis' 
in Paper II. 



a star on the basis of its B z'K colors. The s econd is the 
figure of merit proposed bv IBrammer "eTall (|2008h : 



[(^phot) 



X 



^99 
■'up 



^99 
^lo 



PAz=0.2 



(2) 



This quantity combines the of the fit at the nomi- 
nal redshift, the number of photometric points used in 
the fit, A^fiit, the width of the 99 % confidence inter- 
val, (zyp — zf^), and the fractional probability that the 
redshift lies within ±0.2 of the nominal value, paz=o.2] 
all of these qua ntities are output by EAZY by default. 
IBrammer et al.i (.2008) have shown that a cut of > 2 — 
3 can remove a large fraction of photometric redshift out- 
liers. 

7.3. Photometric Redshifts — Validation 

In AppcndixIXl we describe both the spectroscopic red- 
shift determinations that we have compiled for objects 
in the ECDFS, and show the 2:phot~2spoc agreement for 



individual 



-spec 



samples. For all 'secure' redshift deter- 



minations, the random and systematic photometric red- 
shift error is = 0.036 and med[Az/(l + z)] = -0.025. 
In comparison to spectroscopic redshifts from the K20 
survey, which is highly spectrally complete in the magni- 
tude regime in which we are operating, the random error 
is (Tz = 0.033, with an outlier fraction of less than 5 %. 
(Here, we define the outlier fraction as the relative num- 
ber of sources for which Az/{l + z) > 5az.) We also draw 
particular attention to the excellent agreement between 
our photometric redshifts an d the spectroscopi c deter - 
minations for the sample of IVan der Wei et al.l (|2005| ). 
which is a sample of 28 early type, red sequence galaxies 
at z ~ 1; we find <Jz = 0.022, with no outliers, and essen- 
tially no systematic offset. For comparison, the overall 
photometric redshift accuracy of the COMBO-17 survey 
for our Zspec comparison sample, but limited to Zgpcc < 1, 
is CTz = 0.020. 

However, we also show in Appendix El that none of the 
available spectroscopic samples is particularly represen- 
tative of the MUSYC ECDFS sample. In particular, in 
almost all cases there is a correlation between redshift 
security and (J — K) color, such that redshift determi- 
nations for blue galaxies tend to be more secure, and 
so these galaxies are over-represented among MUSYC 
ECDFS galaxies. Even the K20 sample, which is 92 % 
complete for X^^'^sa) <; 20, does not probe the reddest 
galaxies in our sample, presumably because they are too 
rare to be found in that survey's rather small area. There 
is, therefore, the very real danger that looking only at the 
■^spoc — z-o iiot agreemen t provides a false sense of security 
(see also IBrammer et al.ll200^ . since there are compar- 
atively few ZspecS available for the faintest and reddest 
galaxies in the catalog — especially given that these are 
the main objects of interest. 

For this reason, we have compared our photometric 
redshifts to those froi n COMBO-17 (I Wolf et all [2Ciol 
and GOODS-MUSIC (jCrazian et aLll2006n : the results 
of this comparison are shown in Figure [Ml While these 
comparisons are extremely useful for identifying system- 
atic differences between different Zphot solutions, with- 
out spectroscopic redshifts as a referent, they cannot be 
used to decide which is 'better' in the case of a dis- 
agreement. To this end, the red lines in each panel 
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of this Figure show the spectroscopic redshifts (where 
available) by connecting the (zphot, ^^phot) point to the 
point (^spcc, ^spoc)- In each panel, vertical lines thus in- 
dicate where the COMBO-17 or GOODS-MUSIC Zphot 
is 'right', while the MUSYC Zphot is 'wrong'; conversely, 
horizontal lines show where the MUSYC -Zphot is 'better' 
than that from COMBO-17 or GOODS-MUSIC. Note 
that for the comparison to COMBO-17, we restrict our 
attention to those galaxies with R < 24, since this is the 
rehability Hmit of the COMBO-17 catalog. 

Owing to its medium-band photometry, the COMBO- 
17 redshifts should be significantly better than our own 
for < 1, but without NIR photometry, the redshifts of 
z > 1 galaxies are poorly constrained. The agreement 
between the COMBO-17 photometric redshifts and our 
own (left panel of Figure , the agreement is indeed 
very good for Zphot < 0.8. For i? < 24 and 2:017 < 
1.0, the random scatter between the COMBO-17 and 
MUSYC photometric redshifts is CTz = 0.034; separately, 
for i? < 24 and Zgpcc < 1, the photometric redshift error 
is CTz = 0.030 for MUSYC, and 0.020 for COMBO-17. 

There are, however, several important differences be- 
tween the MUSYC and COMBO-17 redshifts. First, 
note the effect of the zci7 < 1-4 grid used by COMBO- 
17; coupled with their method of assigning redshifts 
(viz., marginalizing over the redshift probability distri- 
bution), this means that galaxies are essentially never 
given ZC17 > 1-3. 

The exceptions to this rule are those objects that 
COMBO-17 has classified as QSOs on the basis of their 
optical SEDs; where MUSYC tends to place these ob- 
jects at 2;mus ^ 1: the COMBO-17 redshifts are very 
good. (Note that we have made no attempt to explicitly 
accommodate AGNs or QSOs in our photometric red- 
shift calcuation.) In the left panel of Fig ure ITU w e mark 
X-r ay selected ga l axies from the iSzokolv et al.l (|2004l ) 
and lTreister et al.) ()2008[ ) catalogs with a cross. For this 
i? < 24 sample. X-ray selected sources are roughly three 
times as likely to be outliers (here, we define outliers as 
those objects with |A2;|/(1 -|- ^spoc) > 0.1): the outlier 
fraction for X-ray sources is 35% (75/217), compared to 
11% (164/1438) overall. Said another way, roughly half 
of all (i? < 24) outliers are X-ray sources. 

Secondly, there are two populations of objects with 
Zen ^ 0.2 that are placed by MUSYC at either zmus ^ 
0.4 or zmus ^ 1-4. From this first population, no ZspccS 
are available; for the second, the handful of available 
ZgpocS confirm that these galaxies are at z > 1.4. On 
the other hand, for the diffuse cloud of galaxies given 
zmvs ^ 0.5 and 0.5 < zci7 < 1.0, the ZgpccS support the 
COMBO-17 determinations. 

Thirdly, while objects given 0.8 < Zphot ^ 1.2 in one 
catalog generally lie in the same redshift interval in the 
other, there is only a very weak correlation between the 
redshifts within this interval: the implication here is that 
objects with Zspec ^ 0.8 are assigned 0.8 < zci7 ^ 1-2 
more or less at random on the basis of optical data alone. 
In other words, while the COMBO-17 Zphot — -^spoc agree- 
ment is excellent for Zgpcc ^ 0.8, a 0.8 < zci7 ^ 1.0 se- 
lected sample may suffer significant contamination from 
Zspec ^ 1 galaxies with poorly constrained redshifts. 

Looking now at the comparison with the GOODS- 
MUSIC redshifts (right panel of Figure [Tl)) . it is clear 



that, while the random scatter between the two deter- 
minations is larger than for the previous comparison, at 
least for z < 1, there are no signs of major systematic 
discrepancies. The random scatter between the GOODS- 
MUSIC and MUSYC photometric redshifts is = 0.065; 
separately, for the same Zgpoc comparison sample, the 
random errors are Cz = 0.036 for MUSYC, and 0.043 for 
GOODS-MUSIC. Both MUSYC and GOODS-MUSIC 
suffer from catastrophic failures, where Zgpcc "^0.7 galax- 
ies are given Zphot ~ 0.2; although this appears to be 
a greater problem for MUSYC. GOODS-MUSIC also 
seems to have some systematic issues for Zphot ~ 0.4. 

In this panel, we mark with a circle those objects with 
Qz > 3. Whereas roughly half (938/1787) of the ob- 
jects plotted in this panel have robust ZgpocS, the frac- 
tion among those with > 3 is just 33% (242/735); 
again, this underscores the importance of having a rep- 
resentative spectroscopic comparison sample. Using the 
cut |(zMUS - ^GDs|/(l + 2GDs) > 0.2 to quantify the 
level of disagreement between the GOODS-MUSIC and 
MUSYC redshifts, objects with Qz > 3 are twice as 
likely to be outliers: the fraction is 60% (99/166) for 
Qz > 3, compared to 33% (586/1787) overall. We note 
that the fraction of sources with Qz > 3 increases from 
< 5% for ZGDS < 1-2 to - 15 % for 1.2 < zqds < 2.2. 
For ZGDS > 2.5, roughly half (9/21) of all galaxies have 
Qz > 3. Similarly, X-ray-selected galaxies are more likely 
to be outliers: the outlier fraction for X-ray sources is 
43% (16/37). 

Again, we caution that, without spectra for a large, 
representative subsample of the objects common to these 
two catalogs, it is not possible to determine whether 
one catalog is truly 'better' than the other. Moreover, 
given the differences between the MUSYC and GOODS- 
MUSIC catalogs — particularly the inclusion of ACS and 
IRAC imaging in the GOODS-MUSIC catalog— it is not 
possible to say whether any differences in photometric 
redshifts are due to the photometric redshift algorithms 
or to differences in the data themselves. Given these 
differences, however, the broad agreement between the 
MUSYC and GOODS-MUSIC ZphotS, and especially for 
•^phot ^ 1 where ZspocS are increasingly hard to come by, 
is certainly encouraging. 

8. INTERPOLATING RESTFRAME PHOTOMETRY 
— INTRODUCING InterRest 

Given an SED and a redshift, we have derived rest- 
frame photo metry following the rn ethod described in Ap- 
pendix C of iRudnick et al.l ()2003[ ). This method is best 
understood as interpolating between two points in the 
observed SED to come up with a restframe flux. We have 
developed an IDL implementation of this algorithm for 
interpolating restframe photometry, dubbed InterRest. 
InterRest has been specifically designed to dovetail with 
EAZY: it accepts the same inputs and configuration files, 
and uses the same algorithms for integration, etc. We 
have made this utility freely available to the astronomi- 
cal community. 

The essential idea is to use a set of template spectra 
to construct a color-color relation for galaxies at a given 
redshift. Specifically, we relate a color in terms of two 
observed filters to another color in terms of an observed 
filter and the desired restframe filter. For example, in 
order to find the restframe r flux of a galaxy at z = 1.2 
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Fig. 15. — Illustrating the InterRest algorithm for interpolating restframe fluxes — Note that under normal circumstances, in order to 
interpolate a restframe flux, we would relate an observed-minus— observed color to a restframe— minus— observed color; in this case we are 
using the {V — I) color to predict the flux in t he observers' R band. By comparing the interpolated and observed R band fluxes, we will then 
be able to validate the algorithm (Figure [TCt . The algorithm works as follows: using a set of template spectra (red crosses), we construct a 
(redshift-dependent) color-color relation for galaxies {red line); once the galaxy color— color relation has been defined, it is possible to read 
off the 'unknown' color (in this case, R — I) of any object, given its known, observed color (in this case, V — I). In both panels, the points 
with error bar s sh ow galaxies in a narrow spectroscopic redshift range, with colors measured to better than 0.05 mag; these galaxies are 
used in Figure [T6l to validate our restframe color determinations. (See main text for further discussion and explanation.) 



(Acm = 6220 A; Aob = 13700 A), we would relate the 
(z' — J) color to the (r^^i 2 — J) color; the r^^i 2 flux 
then immediately follows. 

This process is illustrated in Figure 1151 with one cru- 
cial difference: whereas normally, in order to interpolate 
a restframe flux, we would relate an observed-minus- 
observed color to a restframe-minus-observed color, in 
this example we are concerned with using the observed 
(V — I) color to predict the observed {R — I) color, and 
so the observed R flux. In this way, we will be able to 
test the accuracy of the algorithm, through comparison 
between the predicted and observed R fluxes. Even so, 
the example still serves to illustrate the idea behind the 
algorithm. 

In each panel of Figure [151 the points show the ob- 
served VRI colors of galaxies with spectroscopic red- 
shifts in a narrow interval; we have selected the two most 
prominent redshift spikes, and restrict our attention to 
galaxies with colors measured to better than 0.05 mag. 
The red crosses in each panel show the synthetic VRI col- 
ors for the default EAZY/InterRest template set, which 
we use to construct an approximate color-color relation 
for galaxies at each redshift. In both panels, the default 
EAZY/InterRest template spectra can be seen to do a 
reasonable job of describing the true color-color relation 
for galaxies at each of the two redshifts in question. 

Now, for any individual galaxy (red point, circled; cho- 
sen at random), using the {V — I) color, it is possi- 
ble to read off the (i? — /) color (grey lines) from the 
synthetic color-color relation.. Again, under normal cir- 
cumstances, we would be relating an observed-minus- 
observed color to an restframe-minus-observed color; our 
interest here is in validating the performance of the al- 
gorithm. 

As a single algorithmic detail, it is possible that the 
known-known colors (i.e., {V — I) in the above exam- 
ple) of two templates are very close, but for quite dif- 



ferent known-unknown colors (i.e., {R — I) above): in 
this case, small changes in color or redshift can produce 
very large changes in the final result. To avoid this sit- 
uation, where the known-known colors are too close, we 
simply replace these points with their mean (in magni- 
tude space). This can be seen in Figure fTSl where the 
crosses show the points for the individual template spec- 
tra, and the squares show the points used to construct 
the color-color relation. Algorithmically, we define 'too 
close' as two points being separated by less than 5 % of 
the range spanned by all template spectra. 

In Figure I16[ we show the differences between the R 
fluxes interpolated as described above, and the observed 
R fluxes in the MUSYC catalog, plotted as a function of 
(left to right), spectroscopic redshift, observed color, and 
restframe color. These plots are based on the Zspoc com- 
pilation used in Figure [141 f^^id described in Appendix [Al 
but limited to those galaxies with VRI colors measured 
to better than 0.05 mag. The black points are for indi- 
vidual galaxies; the red error bars show the mean error 
and random scatter in bins. 

Both the random scatter and the systematic offset be- 
tween the observed and interpolated R fluxes are at the 
level of 0.05 mag. There are clear systematics with red- 
shift (i.e. restframe wavelength), which appear to be re- 
lated to the 4000 A break. There also appears to be a 
problem at the level of 0.05 mag for the reddest galaxies 
(u — r > 2). The random error in the interpolated R 
fluxes is typically ^ 0.05 mag. This is comparable to the 
uncertainties in the photometry itself, but probably at 
least partially reflects the intrinsic width of the galaxy 
color-color(redshift) sequence; if so, this represents a fun- 
damental limit on the accuracy of the algorithm. 

Note that whereas we would typically use two neigh- 
boring fllters to interpolate restframe photometry, the 
wavelength span here is roughly twice as large; we there- 
fore expect the true systematic errors in restframe fluxes 
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Fig. 16. — Validating the InterRest algorithm for interpolating rest-fra me photometry — In each panel, we show the difference, AiJ, 
between the R band magnitude inferred from the (V — I) color as in Figure [TSl and that directly observed, plotted as a function of (left to 
right) spectroscopic redshift, observed color, and rest-frame color. The points in each panel show galaxies with robust spectroscopic redshift 
determinations, and VRI colors measured to better than 0.05 mag; the red points with error bars show the biweight mean and scatter 
in AR. The random error in the interpolated R band magnitude is typically < 0.05 mag; comparable to the observational uncertainties 
themselves. Systematic uncertainties, as functions of both redshift (i.e. restframe wavelength) and restframe color are at the level of < 0.05 
mag. Note, however, that the (logarithmic) wavelength interval between V and I is roughly twice as large as we would normally use to 
derive rest-frame photometry for real galaxies. We therefore present these numbers as upper limits on the true errors; we expect the true 
errors to be smaller by a factor of 2 — 4. 



TABLE 5 

Summary of the Contents of the Restframe Photometry Catalog 



Column No. Column Title Description 



1 




id 


Object identifier, beginning from 1, as in the photometric catalog 


2 




redshift 


Assumed redshift; we use either the z_m2 value output by EAZY, or the spectroscopic red- 
shift, where available. 


3— 


-17 


RF_F1, etc. 


Restframe photometry for Bessel U BVRI filters" 


18- 


-32 


RF_F6, etc. 


Restframe photometry for Johnson-Cousins UBVRI filters 


33- 


-47 


RF_F11, etc. 


Restframe photometry for Gunn ugriz filters 


48- 


-54 


RF_F16, etc. 


Restframe photometry for GALEX NUV and FUV filters 


55 




distmod 


The distance modulus implied by redshift, assuming a given cosmology'' 



Note. — ° For each object and filter, InterRest outputs two fiags; extrapra (where n refers to the restframe filter number), which 
indicates where it has extrapolated beyond the observed SED, and widegapn, which indicates where it has not used neighboring filters due 
to, for example, missing or negative photometry. ^ Note that the fluxes output by InterRest are observed fluxes through restframe filters; 
that is, they have the same units as the observed, input photometry. The user must therefore perform the conversion to apparent and 
restframe magnitudes using the appropriate zeropoint and distance modulus. 



(c/. colors) to be 2 — 4 times smaller than in the above 
example; i.e. at the level of 0.01 — 0.02 mag. 

As a final aside, we note that we acheive comparable 
accuracies using the E, Scd , Sbc, and Im templates from 
IColeman. Wu fc WeedmanI (Il980l), su pplemented with a 
starburst template from iKinnev et al. (1996.) • These 
templates are plotted in Figure [15] (black squares, la- 
beled) for comparison to the defau lt EAZY/InterRest 
templ ates. We have also tried using iBruzual fc CharlotI 
([2003') synthetic spectra, assuming Single Stellar Popu- 
lations (SSPs; logt = 6.5, 7.0, 10.0, 10.3 Gyr) with 
a Salpeter IMF and solar metallicity, and no dust ex- 
tinction (shown by the dashed yellow line in Figure [T5|) . 
Using BC03 spectra, we find serious systematic errors — 
on the level of up to 0.2 mag — both as a function of 
redshift, and of restframe color; this is true whether we 
assume a SSP or exponentially declining star formation 
history. These models do not reproduce the observed 



colors of real galaxies, an d so are unsuitable for this pur- 
pose. Similarly, using the lBlanton fc RoweisI (|2007i) tem- 
plate set, which are derived from a library of BC03 spec- 
tra with a wide range of ages and metallicities using the 
NMF algorithm, we find peak-to-peak systematic errors 
at the ~ 0.1 mag level; the random errors are also at the 
0.1 mag level. 

In Table O we summarize the contents of the restframe 
photometry catalogs that we are releasing: note that we 
provide two separate catalogs based on photometric and 
spectroscopic redshift determinations, respectively. 

9. SUMMARY 

We have described a new i^T-selected catalog of the 
ECDFS based on existing optical and NIR data, supple- 
mented by original z'JK imaging taken as part of the 
MUSYC project. The final UUssBVRIz' JHK photo- 
metric catalog (gl Table [2|) covers 900 □ " to a (5(t, 
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point source) limiting magnitude oi K = 22.0 mag; note, 
however, that H band data is available for only 80 % 
of the field. Included in the photometric catalog are a 
spectroscopic redshifts for 2914 unique objects, collected 
from the literature (Appendix E]) . In addition, we are 
also making available a photometric redshift catalog, de- 
rived from the MUSYC ECDFS photometry using EAZY 
( H7.21 Table S]), as well as catalogs of interpolated rest- 
frame photometry generated using InterRest (fJSl Table 

The data described in this paper will form an impor- 
tant part of two ongoing NIR survey projects. The K 
imaging is key for analysing the SIMPLE IRAC data 
(Damen et al., in prep.). The broadband imaging pro- 
vides the backbone for an optical medium-band survey, 
which will add 18 additional bands (Cardamone et al., 
in prep.). There is also a NEWFIRM medium band NIR 
survey planned, which will allo w much greater photomet - 
ric redshift accuracy for 2 > 1 ()Van Dokkum et al.|[2009l ). 
We have invested significant time and effort in validating 
the absolute and relative calibration of the imaging data, 
as well as our analysis techniques, so as to maximize the 
legacy value of our catalogs. We summarize the results 
of these checks below. 

Astrometry — The relative astrometric calibration of 
each band has been validated to O'.'IS (0.56 pix). In ab- 
solute terms, the absolute astrometry is accurate to O'.'S 
(1.12 pix), with a slight shear across the field at the level 
of O'.'l (0.37 pix; see §??). 

Completeness — We have quantified the completeness of 
the catalog for sources with an R^/^ profile in Figure HJ 
we present these values as lower limits on the complete- 
ness. While the catalog is formally surface brightness 
limited, a comparison to much deeper NIR imaging over 
the GOODS area of the field suggests that the catalog is 
more nearly flux limited. This comparison suggests that 
for K = 22, the catalog is ~ 85 — 90 complete, and > 95 
% rehable (g21). 

Photometric Calibration — While there are significant 
differences between the photometry in the COMBO- 17 
and MUSYC catalogs of the ECDFS ( ^5X2]) . a compar- 
ison between the MUSYC and GOODS photometry in 
the region of overlap validates the MUSYC photometry 
to < 0.05 mag fi )5.2.1[) . We have refined the basic pho- 
tometric calibration using the observed SEDs of main 
sequence stars; we estimate that after this recalibration, 
the photometric cross-calibration is accurate to < 0.02 
mag fi J5X3l) . 

Photometry — Random and systematic photometric er- 
rors due to various aperture effects (including astrometric 
errors and imperfect PSF matching) are limited to < 0.03 
mag and < 0.006 mag, respectively (Figure [3|). We have 
applied corrections to SExtractor's AUTO flux measure- 
ments to account for missed flux and background over- 
subtraction; for synthetic i?^/^-law sources, these correc- 
tions typically reduce the offset between the known and 
recovered total fluxes by 0.05 — 0.10 mag ( §4.3p . We have 
also demonstrated that the photometric errors given in 
the catalog accurately trace variations in the background 
RMS in the NIR images (Figure E]). 

Spectroscopic Redshifts — We have collected and col- 
lated 5374 spectroscopic redshift determinations from 



literature sources, of which 3815 are matched to 2914 
unique sources in our catalog (Appendix |A|). Of these, 
2213 redshifts are deigned 'secure', including 247 stars, 
and 1966 z ^ galaxies. 

Photometric Redshifts — There are some systematic dis- 
crepancies between the COMBO-17 and MUSYC pho- 
tometric redshift determinations in the ECDFS, ow- 
ing to the lack of NIR data in the COMBO-17 cata- 
log; where available, spectroscopic redshifts validate the 
MUSYC values. The agreement between the MUSYC 
and GOODS-MUSIC photometric redshifts is very good, 
however there are a significant number of catastrophic 
errors in both redshift catalogs (Figure [T4|) . In com- 

J iarison to spec trosco pic redshifts fr om the K20 survey 
Cimatti et al] l2002t iMignoh et all |2005^, the random 
photometric redshift error is CTz = 0.033, with an out- 
lier fraction of 4.7 %; the outlier fraction is significantly 
higher for X-ray-selected spectroscopic redshift catalogs 
(Appendix El . 

Restframe Colors — We have interpolated restframe 
photometry for the galaxies in our catalog using an IDL 
utility called InterRest (Sj8]); we also make this utility 
publicly available. Estimated systematic errors in these 
interpolated restframe fluxes, as functions both of rest- 
frame wavelength and of galaxy color, are estimated to 
be < 0.02 mag (Figure fT6|l . Random errors inherent to 
the algorithm are at a similar level. 

The primary science application of the if-selected 
catalog that we have presented here is to characterise 
the properties of massive galaxies at z < 2, including 
their evolution. In Paper II, we demonstrate that this 
catalog is approximately complete (volume limited) for 
M^, > 10^^ M0 and Zphot ^ 1-8, and use this catalog 
to quantify the z < 2 evolution in number density and 
color of massive galaxies in general, and of red sequence 
galaxies in particular. 

In this context, the MUSYC ECDFS dataset provides 
a valuable complement to existing optical surveys in 
the ECDFS targ eting the z < 1; e.g. the COMBO-1 7 
(|Wolf et al.ll200l and GEMS projects (iRix et al.ll200l . 
Further, the z < 2 comoving volume contained with 
the ECDFS field is approximately three times greater 
than that at z < 3.5 within the GOODS region in the 
CDFS. The MUSYC ECDFS catalog thus also comple- 
ments the much deeper GOODS-CDFS data, by allow- 
ing better sampling of rare objects, including the most 
massive galaxies at moderate- to high-redshifts. Taken 
together, these combined datasets form an outstanding 
laboratory to study the basic properties of galaxies over 
nearly 90 % of the history of the universe. 
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APPENDIX 

A. A COMPILATION OF PUBLIC SPECTROSCOPIC REDSHIFT DETERMINATIONS FOR THE 

MUSYC ECDFS CATALOG 

The ECDFS has been targeted by a number of large spect roscopic redshift campaigns , including : optical spectroscop y 
of the origi nal CDFS X-ray catal og bv lSzokolv et al.l (I2004D. the K20 survey (Cimatti eFaI] |2002l: iMignoh et aIll2005D. 
the V VDS (|Le Fevre et al.l[200l . the GOODS F0RS2 CVanzejlaelaLI [20051, l2006i, i2007D and VIMOS (IPopesso et all 



[200i) ca mpaigns, the IMAG ES survey (jRavikumar et a l. 2007) , a MUSYC program targeting X-ray sources in the full 
ECDFS (|Treister et al.l[200l . and a VIMOS campaign by S K oposov et al.l (|in prep] ). A summary of the spectroscopic 
redshift resources we have used is given in Table El Altogether, we have collected 5374 separate spectroscopic redshift 
determinations, of which 3815 are matched to 2914 unique objects in our catalog. 

In cases where multiple spectroscopic redshift determinations/identifications are available for individual objects, our 
guiding principles for selecting a redshift were as follows. First, we adopt the most common redshift determination 
(where Az < 0.01 is taken as agreement, and we do not consider repeat observations by the same team as an 
independent measurement). 574 objects in the catalog have multiple, consistent redshift determinations. Where 
there is no consensus, we discriminate between redshift solutions on the basis of the figure of merit developed 
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TABLE 6 

Summary of the Spectroscopic Redshifts Available for MUSYC ECDFS detections 



Reference(s) 


Source 


Internal 


No. 


No. 


Median 


A T TV T A 

NMAD 


Outlier 




Code 


Qual. Flag 


Galaxies 


Adopted 


Az/{1 + z) 


Az/{1 + z) 


Fraction 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Cimatti ct al. (2002); 


K20 


1 


267 


232 


-0.025 


0.033 


0.047 


^gnoli ct al. ('2005') 







14 


2 


-0.012 


0.069 


0.182 


Szokolv et al. ('2004') 


Xray 


> 2.0 


114 


114 


-0.024 


0.037 


0.135 






< 2.0 


17 


4 


0.045 


0.146 


0.133 


Le Fevre et al. (2004') 




4 


172 


131 


-0.030 


0.027 


0.027 






3 


347 


267 


-0.030 


0.032 


0.035 






2 


342 


19 


-0.022 


0.058 


0.080 






1 


82 


1 


-0.003 


0.127 


0.017 






9 


49 


1 


0.016 


0.199 


0.036 


■Vanzella et al. (2005. 2006. 2007') 


GDS-F 


A 


306 


226 


-0.023 


0.044 


0.034 




B 


77 


14 


-0.029 


0.080 


0.054 






C 


52 


4 


0.025 


0.106 


0.079 


Pooesso et al. (2008") 


GDS-V 


A 


289 


197 


-0.036 


0.030 


0.048 




B 


59 


3 


-0.026 


0.081 


0.087 






G 


48 


1 


-0.008 


0.144 


0.051 


Ravikumar et al. (2007") 


IMAGES 


1 


267 


219 


-0.032 


0.030 


0.067 






2 


168 


24 


-0.025 


0.046 


0.056 






3 


51 


7 


-0.012 


0.095 


0.000 


ireister et al. (2008) 


MUS-I 


N/A 


165 


120 


0.001 


0.112 


0.125 




MUS-V 


N/A 


34 


33 


0.011 


0.295 


0.000 


S KoDosov et al. (in Drep.) 


Kopsv 


N/A 


455 


283 


-0.034 


0.025 


0.043 


Groom et al. (2001") 


KX 


N/A 


17 


5 


-0.016 


0.029 


0.353 


Strolger et al. (2004') 


SNe 


N/A 


9 


2 








Van der 'Wei et al. (2004. 20051 


vdWel 


N/A 


28 


26 


-0.007 


0.022 


0.000 


Daddi et al. (2005') 


Daddi 


N/A 


5 


5 








Dohcrty ct al. (2005) 


LCIRS 


1—3 


14 


10 


0.003 


0.050 


0.071 


Krick ct al. (2006) 


Krick 


N/A 


12 


12 


0.056 


0.134 


0.000 


Total 






2863 


1966 


-0.029 


0.036 


0.078 



Note. — For each spectroscopic rcdshift sample we have used, we give both the redshift source catalog (1) and the identifier used in the 
MUSYG Zspcc catalog (2); further, we have broken up each sample by the internal quality flag (3), where available. For each (sub)sample, 
we give the number of galaxies matched to the MUSYC EGDFS catalog (4), and the number of galaxy redshifts adopted in the final catalog 
(5). We also give the systematic (6) and random (7) photometric redshift error, computed as the median and NMAD of A2/(l + z), and 
the outlier fraction (8), defined as the fraction of galaxies with Az/(1 -|- z) > 0.1; these quantities are all computed for galaxies in our main 
scientific sample (i.e. those galaxies counted in column 4 with coverage in optical and NIR bands, and with X < 22 and K S:N > 5). 

by iBrammer et al] (|2008f ). evaluated for the spectroscopic redshift. An exception to this rule is for redshifts from 
the X-ray selected catalogs, which do occasionally have extremely high values of Qz{zs-pcc)i even when confirmed by 
other secure determinations from other catalogs. Where Qz{zspec) does not clearly discriminate between the possible 
solutions, we fall back onto the quality flags given by the different spectroscopic surveys. Note that for this purpose, 
we do not consider the VVDS '2' flag as 'secure'. Similarly, we give preference to the results of smaller studies, which 
presumably have devoted greater care on a per object basis. Reassuringly, in almost all cases, these criteria reinforce 
one another. Finally, we choose to adopt redshifts from sources that provide classification information where available; 
this means that we tend not to adopt redshifts from, for example, the VVDS catalog where other determinations are 
available. Moreover, we consider X-ray selection as an additional piece of classification information; accordingly, we 
adopt redshifts from the Szokoly et al. (2004) and Trcistcr ct al. (2008) catalogs where available. 

In this way, we have constructed a compendium of spectroscopic redshift determinations for 2914 unique objects in 
the MUSYC ECDFS catalog, including 283 spectrally-classified stars. Although all of these determinations are given 
in the catalog, we will only consider those deigned 'secure', either by virtue of their quality flags, or through agreement 
between multiple sources. This leaves 2213 robust spectroscopic redshifts for objects in the MUSYC catalog; 1966 of 
these objects are identified as ^ ^ galaxies. 

Figure [17] shows the Zphot — ^spcc diagram, broken up by the Zgpoc source catalog, and quality flag. Within each 
panel, we give the NMAD and median offset in Az/(1 -I- z)] these values are also given in Table El the grey region 
indicates the Sctz errors around the Zphot — -^spcc line. Above each Zphot — -^spcc diagram, we also show the distribution 
of each Zgpoc sample in observed (J — K) — K color — magnitude space, in comparison to the full MUSYC catalog. 

For 'secure' redshift determinations, the 2;phot~-2^spec agreement is really quite good: the typical random scatter is 
(Tz < 0.040. Particularly for Zgpoc ^ 1, we do appear to slightly underestimate galaxies' redshifts; typical systematic 
errors are Az/(1 + z) ^ —0.025. For the iSzokolv et all (j2004) catalog, the random scatter in Zphot determinations 
is still quite good, but for the MUSYC spectroscopic redshift program (|Treister et al.l l2"008f ) . which targets brighter 
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Fig. 17. — The Zphot — ^spec diagram for individual spectroscopic redshift sources and quality flags — We show the agreement between 
the MUSYC photometric redshift and many literature spectroscopic redshift catalogs, broken up by internal quality flag where available. 
Within each panel, we give the random photometric redshift error, Cz, as well as the outlier fraction, which is defined as the fraction of 
objects with Az/{l + z) > 5crz. The grey regions in each panel indicate the 3(Tz regions. Above each Zphot — Zspcc diagram, we highlight the 
particular 2:spoc sample in (J — K) — K color— magnitude space. It is harder to obtain robust Zspec determinations for redder galaxies; these 
galaxies are therefore underrepresented in all Zspec samples. For this reason, we have validated our photometric redshift determinations 
through comparison with those from COMBO-17 and GOODS-MUSIC (see Figure [Tit. 



X-ray sources, the Zphot — Zgpcc agreement is poor. 

Further, while the outher fraction is generaUy at the level of a few percent, catastrophic redshift failures appear 
to be a significant problem for X-ray selected sources, (Recall that we make no attempt to explicitly incorporate 
AGNs or QSOs in our photometric redshift calculation,) Among X-ray-selected sources, the fraction of galaxies with 
|Az|/(l-|-z) > 0,15 is 30% (82/271); for the full gamut of robust spectroscopic redshifts, the fraction is 9% (178/1966), 
Said another way, 46 % (82/178) of all outliers are X-ray sources. 

We draw particular attention to the comparison with the results from K20, which is highly spectrally complete in 
the magnitude range that we are operating in. In comparison to the K2 redshifts, we hav e achie ved a photometric 
redshift accuracy of — 0,034. We also draw attention to the sample of IVan der Wei et al,l ()2005| ). which consists of 
28 early type, red sequence galaxies at z ~ 1, for which we have achieved a photometric redshift accuracy of az = 0,022; 
in fact, this is the sample for which we have the best photometric redshift agreement. 

The crucial point to be made from Figure 1171 however, is that since most of the different Zgpcc samples that are 
available in the ECDFS are not NIR~selected, they are not generally representative of the sources in our photometric 
catalog. For this reason, we validate our photometric redshift determinations in H7,3I through comparison with the 
COMBO-17 and GOODS-MUSIC photometric redshifts. 



